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CAVITATION DAMAGE OF ROUGHENED CONCRETE SURFACES* 


Donald Colgate! M, ASCE 


he ‘paper discusses an exploratory iensettaire study concerning the eval- e 


uation « of the cavitation potential of various types of roughened concrete sur-— 
: faces. _ Methods of evaluating the data determined by the laboratory study, veil 


vigorous ‘campaign cavitation erosion of hydraulic 
_ continues } Space. . Technical personnel in n the > rand laboratories of o our uni- = 


oo and unrelated plausible explanations as to just how this harmless 
appearing vapor | cloud can inflict such unbelievable damage to construction — i 
materials. The subsequent discussions that are encouraged when these aiffer- 
nt ities are aired before all who are interested in the problem have led to | 


3 * he remarkable advances made thus far in the methods of protecting « our hy 
‘The common construction material 1 most readily damaged 


e varied. However, due to its susceptibility to cavitation erosion, the pn : 
tions configuration of, and the very texture of, concrete 


oss the « concrete to “hump” about. 3 inches in at least one of the 5- foot 


lifts. An examination of the spillway face after about 9 seasons of operation 
a 


= emai open until April 1, 1960. To extend the closing date one month, a 

_ written request must be filed with the Executive Secretary, ASCE. Paper 2241 is 

part of the copyrighted Journal of the Hydraulics Division, Proceedings of th - ae : 
_ American § Society of Civil ‘Engineers, Vol. 11, No. HY 11, November, 1959. 
a. Presented at the October, 1956, ASCE Convention in Pittsburgh, a 
Engr., Bureau of Reclamation, Denver Federal 
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disclosed that tl the concrete downstream trom. the hum ) had been dam 


= cavitation (Fig. 1). In this instance the « designers had specified that the spill- 
a way section conform in all respects to the “lines, grades, and dimensions” 


4 expected, but the permissible limits of such ‘deviation to prevent 
_ Cavitation have not been accurately determined, therefore, in light of the __ 
3 Grand Coulee experiences, extremely rigid specifications are now adhered to. 


eh It is common practice to evaluate extensive damage where cavitation is we 


the large eroded areas downstream aret the result of jet action. From a prac 
Sem viewpoint, the ‘mechanics of the total damage is of secondary et 


_ The relative destructiveness of a he and of ‘cavitation has ‘been inadvert- 


Fig. 1 Grand Costes: Shov 
and Cavitation Erosion. 
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sible, surfaces an avitation pocket. This does 
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— remain and set up; merely sti 
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‘CAVITATION DA DAMAGE 


“not these roughened surfaces would induce cavitation. Plaster molds were a. 

_ made of several different damaged areas, and concrete casts made from these 
Lah os ‘molds, so that the cavitation potential of the various surfaces might be cain 
2 ated in the Bureau of Reclamation Hydraulic Laboratory. However, spillway | 
lenses were imminent and insufficient time remained for a laboratory study, 

Ps therefore, the repair of all the scoured area to the original contours was 
li ae The problem concerning the severity of surface roughness which might be 
ae tolerated remained to be solved. Since the casts of the roughened surfaces “a 


1 7 from Davis Dam were on hand, the Bureau initiated a limited exploratory pro- a 
gram to devise a means of evaluating various roughened surfaces as regards 


their cavitation potential. program envisioned eventual classification of 
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Fig. 2: -foot- -square Plaster Mold of Eroded Surface of I 


| about until quite an area was scoure 
ore ini some places to a depth o q 
pre below the finished surface (Figs. 2 and 3). When: depth of Sinches or 
inished surface (Figs. 2 and 3). When the bucket was 


surface on new installations for r repairing 
he laboratory testing apparatus was designed and built to accommodate the ~ 


Davis Dam specimen. apparatus (Fig. 4) permitted a stream 6 inches| 


was made of a transparent material so that the presence of cavitation could ae are a 
? _ be determined visually (Fig. 5). During the preliminary test run it was a 
termined that, under certain conditions, cavitation did form on the protruding _ 


_ and that of Specimen No. 3 about 1/4 in 


the discharge piping was so arranged that the pressure head could be lowered 
34 to about 17 feet of water below atmospheric pressure . The measurements — 


a a In the analysis the data are applicable to either closed condui “ a 
_ channels. Although the laboratory study was made in a closed conduit rectan-— 
Pa ome” in cross section, the computations were ‘made using circular pipe i 
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=a cimen No. 3). 
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co 
fa 
on 
= 
> = F No.3. The average exposed aggregate of specimen | ed 
men No. 2 and Specimen lowest point of the roughened surface, 
out 85 fos through the test section. = 
scharge ressures al Various points On the viewing wh 


shear, and the stream tne: affected by a given s sur- 
eae i = face in the same manner for either closed conduit or open channel flow The 
Ais be computations regarding the velocity distribution in the stream above the ‘sur- Bo 
face recognize that the Prandtl universal logarithmic velocity- 
| pt law for pipes applies equally well for an infinitely wide open channel. -_ 
The analysis ¢ considered that the tested surface was a full scale 
a and that the boundary shear, or shear velocity, of the surface could be deter- "re ‘ 
/ mined directly from measurements made in the laboratory apparatus. A plot a 


on semilog paper of the profile perpendicular to the specimen 


the 


K = a constant related to the boundary 
she: 
) and solving si- 


values. two points (Vi 1> ) and (Wa, Yo 
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or the 


= profile measured above the specimen. . This procedure w was followed _ 
several combinations of velocity and pressure. A plot was then made > ro 
relationship betweeen the shear velocity and the pressure on the viewing 
3 This curve is readily obtainable from the test apparatus, but the ontethti- 
ity is remote that flowing water in a field installation will have been in con- — f 
tact with the roughened surface long enough to establish uniform turbulent — = 
- flow, and an accurate determination of the shear velocity would be virtually 
ee | study of various } velocity profiles (Fig. 9) reveals that, for the same . 


— smooth surface than for a rough one. If the approach to the roughened area 
was reasonably smooth surface sufficiently long to establish uniform flow, 
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v= 40 fos Ave model v= 80 fps 
Pressure heod= Pressure head=+143.| 
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Computes. velocity ‘for smooth ant rough 
in turbulent flow. 


2 
4 gh surface wou e velocity hich mu hened area. 
suriac ince the one w rougnened & tex 
the ro ig. 7). Sin it is the ed or 1 imen No. 
the roug would exist on hened s 
attack vitation hat a roug t = 
her cavitation me tha ith wate — 
= 


this ‘surface | for a distance From 
i the shear velocity at which cavitation would occur would be 4.2 fps. s Consider 
re. the velocity about 0.3 inch from the mean | surface to be critical. ‘This veloc- 
ity may be (for this surface) from the relation: ina 


we from | the average of the 


Vo y 0.3 the mean 
i105 + 27) 594: 


= (0.37)(5) = 1.85 feet from the bottom 

Channels, Civil Engineering, June 1941, 
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he cavitation will occur on Specimen No. 2 
with a relatively smooth approach sur- 
a In th the case of flow with a blunt velocity profile, the average | velocity can r 
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ssmooth- 
ar of the roughened area. 
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cipient cavitation to exist on the roughened | surface in the illustration: — 


No. 2 texture A fps 
Blunt velocity profile 


velocity near nein the labors surface versus depth of flow pressure 
: Py. head) | for the condition of incipient cavitation for the two tested surfaces (Fig. Fs 
ee * A family of curves of this type f for various surface textures would be a % 


xtremely valuable to the designer, or to the field engineer, since the values _ * 
for the chart can be readily computed from known flow conditions cilia wt te 


any surface. There isa need for | the collection, evaluation, | and classifica- 


‘enable the designer } or field engineer - to state with confidence that a a given sur- 

ace texture will or will not induce cavitation 

40 


“TCurves reflect the combination _ 
| of velocity and head producing 
incipient cavitation. Safe | 
operating range is above 
the ‘line. 


Davis Dam.----~ 
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Figure lo. elationsiip 
and average stream velocity, feet per second. 
ae *Fig. 10 shows the most rigorous conditious: that of a stream flowing with a 


blunt velocity profile. = a stream has a ssemalty or fully developed boundar 
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N CHANNEL TREATMENT IN Los. ANGELES “a 


william R. Ferrell, Jr. AM, M. ASCE 


constructed to the f financial burden of flood and for 


untain channe worse in Los Angeles County are undertaken — 
* accomplish three main objectives. These are: (1) to protect the inhabi- 
tants and property below the watersheds from damage by flood - -borne rock 


and debris, (2) to lessen the financial burden of maintaining the large casts. de 
debris - -retaining and flood - -regulating structures that have been constructed 


The population of the Los Angeles Coastal Plain has increased from 
: 3, 800,000 to 5,600,000 persons in the past 10 years and is projected to reach 
a 500,000 by 1970, Keeping pace with this growth has been the increase in ae * | 
Ee housing and industrial development. z The impression 01 one > gets | from an aerial 
view of the broad alluvial cones that have been built. up during ‘the past ae 
thousands of years immediately below the watersheds forming the perimeter 
3 of be Coastal Plain is that here is where most of tl this growth is taking pare. ee. 


4 
“: The San Gabriel Mountains, which form the northern boundary ¢ of this peri- : = 


a meter, extend from the San Fernando Valley in the west, eastward for 55 a 
miles ‘to and beyond the L Los Angeles -San Bernardino County lin line (see location 
map, Fig. 1). With the exception of a few watersheds in the nearby Verdugo Pe 


_ Mountains, watershed treatment efforts have been and will continue to be “aga a 
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the some 116 watersheds of this range which empty 
_ onto the Coastal — watersheds range in size from 0,1 to — 5 
‘The rapid increase in the number of debris basins and dams required to 3 

_ provide protection to this expanding populace from debris-laden flood flows 
has resulted in the impounding an and subsequent ‘removal and disposal of the = 
bulk of the debris that debouches from the watersheds. There are presently | 
= debris basins and debris dams along the San Gabriel front, and this num- 
ber will be increased to 56 by 1970. The over-all flood control system also ~ 
includes 17 flood-regulating dams which are continuously losing valuable 
capacity to impounded debris. Some of these have already been re- 
duced in n capacity to such an extent that mechanical debris removal has been 
: undertaken, In a study(1) recently completed, it was estimated that over oe 
300 million cubic yards of debris could be produced from the frontal iene: 
_ sheds of the San Gabriel range over a fifty-year period of storm activity Wie ee 
. all storms with a frequency of once-a-year to once-in-fifty-years occurred 


| 


during the period . This, Of course, is an extreme 1e assumption but may | be ee 


t 
on. +The situation is made more critical from a financial standpoint in that the — Al 
a ‘cost of handling each cubic yard o of ¢ debris has risen from $0. 33 in 1938 to i 
50 per ‘cubic yard today, and may approach $3.00 per cubic yard by 1970. 
_ These costs, when applied to the estimated debris potential stated above, fi 
= dollar amounts in excess of the financial ee of the Flood Co 


servation of the water supply that is available the that 
: on these watersheds. The value of ground water in this area has already ogo 
‘Feached the vicinity of $20.00 per acre -foot and is expected to approach 


- $45.00 within the next ten years. It should be emphasized that over 50 os all 
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debris being produced from these watersheds, several important = 

projects utilizing check dam systems have been ‘undertaken. | Treatment ef- 

forts have | been based on the theory that if the channels of a watershed can be 

es ‘stabilized, further  downcutting into the e mountain mass will be halted, and the 

major or source of erosion will have been removed. ‘Bank sloughing, which fol- 
_ lows the creation of unstable side slopes due either to downcutting of the chan - a 
nel bottom or to undercutting of f the channel ba banks, is stopped. . The importance a 
of preventing bank s sloughing should be stressed, for | once it starts, it gr adual- B 
_ ly works farther and farther up the hillside contributing large quantities of b q 

. as it progresses, until a stable slope has again been established. In 

many cases, the erosion scar may continue all the way to the ridge line. 
Channel stabilization is also an effective means of halting erosion from 
previously- created slide areas. By providing a Stable toe at an elevated posi- 
: tion upon which the material raveling from above may accumulate, and by re- 
a ducing the erosive power of the stream to wash away this material healing . 
areas may be accomplished. reduction in stream energy occurs as 4 
the result of a 25 to 30 percent ‘decrease in channel gradient caused by the a 
deposition on a milder slope of material behind the stabilization structures. 

5 ‘Further energy and | erosive power is lost through the spreading of the flow 
Es over the widened channel. Once the eroded material has begun to build up at . 

fi by the base of the slide area, bee pean may re-establish itself to further aid 

a Other benefits of channel stabilization are realized through the new debris — ie 

ae and water storage capacity provided by the system. The pore storage in the ee 

_ debris wedge behind each dam is in essence a small reservoir which p> to ‘a 
: _ water conservation in two ways. First, water is stored at time of high runoff aa 
ba: and released slowly at times of low flow, thus correspondingly reducing waste 

_ to the ocean at times of peak discharge. Second, the saturated debris wedge 

- greatly increases the area of water contact in the canyon bottom thereby in- _ 

creasing deep percolation opportunity. Though the salvage in each of these 
reservoirs is minute in itself, the aggregate behind a system of check dams 
in a given canyon is significant. This is a particularly important item in the - i 
_ §San Gabriel Mountains because the ground water storage available in the ere 
- = cracks and crevices of this pe th -fractured — is far | greater than oe Bt 


first stabilization preted of importance constructed in Brand q 
Canyon, near Glendale, in 1938 (see Fig. ‘This project has | proven 
valuable in evaluating benefits from this type of treatment. Nineteendams, 
‘ranging in height from 10 to 40 feet above the stream bed, were constructed om 

. a) stabilize 4,800 feet of channel, Since construction, this system has stored 

_ over 80,000 cubic yards of debris, _ The impounded debris has reduced the 

a average channel gradient through the stabilized reach from an average — 

_ 12.5 percent to 7.3 percent, and has elevated the stream bed from its ian- sg 


 trenched position an average of 11.4 feet through the reach. A typical reach 

ee the Brand Canyon Channel is shown in Fig. 2 along with the plot of i, ; 

of the channel bottom, The profiles were taken | in 1931, 1934, and 


; g 1931 ond 1934 is ‘Clearly shown by the e respective profiles. In the reach below vs 
18+10, however, the stabilizing i influence of two small rock w 
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check dams in 1931 is 1987 the ag ag- 


Poa gradation of the channel resulting from the installation of the large check — 


) a believed that this system has now reached its maximum storage 


capacity. ‘Therefore, a study(2) of debris production for the 1957- -58 storm 
2 season was made to measure | the effectiveness of the system. By coincidence, — 
| _ this proved to be an excellent season to make the study, in that precipitation 
a was nearly 165 ) percent of normal, Sunset Canyon, which is an adjacent water- _ ; 
but with no stabilization treatment, was used for comparison, The study 
was greatly facilitated by the fact that both watersheds have debris- -retaining <a 
= to enable measurement of debris production. Also, physiographic cuales 
features such as elevation, geology, and compass orientation of both water= 
sheds are very nearly the same. Rain gages were located in both ell 
to reflect any significant difference in precipitation a1 amounts. . The drainage a, 
area above Sunset Debris Dam is 0.44 square miles, and that above Brand jot | - 
‘Debris Basin is 1.03.equare miles. 
ee results of this study show that a rate of 8,410 cubic yards per square = 
- mile were eroded from Sunset, while the rate from. Brand was only 50 a 
_ ‘yards per square mile. An analysis of the debris produced from each of the 
_ watersheds indicated that, for Brand, the maximum particle size was 3/8"* in ee 
= diameter; in Sunset, in addition to the finer material, there was a high per- > 
centage of larger material ranging in size up to 1.5 feet in diameter. A study 
_ was also made of erosion source areas in and adjacent to the main channels a i gg 
# of each of the two watersheds. The results showed that in the 3,570 feetof §§— 
_ channel mapped in Brand, 1,850 square feet of active erosion area existed. bases oe 
4 In Sunset, where 3,030 feet ‘of channel were mapped, there were 5,700 square ae 
= feet of ‘active. erosion nm area. ' These two studies indicate that erosion in Brand 
Bs was limited to sheet erosion; in Sunset, in addition to expected sheet we 


-£ the main | Source of debris was the channel and bank areas. oe Thus, it is — 


e tively. The acer of stabilization that has been achieved is shown in Fig. 3. _ i 
by two comparative photographs of the channel and adjacent side slopes Sa 
Be the vicinity of Dam 19 (see Fig. 2). Photo (a) was taken in 1939(3) and ehantt . 

be ‘The e next tr treatment project of interest was completed in in 1948 in the Arroyo 
if Seco watershed, near Pasadena. This project was also experimental in ngture 


Fig D Dam 19 Glendale, Califor nia 


Both photos taken from same point, look- 
were 
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_ reinforced concrete, concrete crib, soil cement, steel bin, and steel arch. gas? ee 
Accurate cost records wi were maintained on all ‘phases of the project to. enable 
- an evaluation of the cost of each type. Although this project provided much Se Gai 
valuable data relative to structure size, design, and desirable construction — 
materials, it proved to be very costly . AS a a result, a damper was placed on — 4 
upstream engineering in this area for a number of years. 
_ By 1956, for reasons previously stated, it became evident that upstream 
engineering was an essential phase of the over-all flood control effort. It was 
therefore decided that a new project should be undertaken that would incor-— is 


porate the experience gained from previous projects i in this area, as well as 


available from similar works in Europe and elsewhere. Following 


thorough study, it was decided that the design and construction of this — 
"should be based on the following general criteria: (1) the spacing of struc- 
ms tures should be such that the impounded debris behind each extends to the toe a 
of the next structure upstream, (2) in general, concrete crib-type construc- ~ 
e tion should be used, (3) the spillway height should not be less than 10 feet —/ 
above the stream bed nor greater than 17 feet, (4) a stabilized debris slope of ea 
0. 7 of natural gradient should be used for the spacing of structures, (5) the | ie 
j structures should be designed to meet all requirements of a gravity-type 
4 “overflow dam, (6) stabilization treatment should be limited to channel gradi- _ a 
_ ents of 20 percent and less, and (7) production and assembly line methods of = 2 
construction should be employed wherever possible. 
Canyon, near La Crescenta, was selected as the watershed to 


‘Control District and the U. S. Forest Service and was completed in late 1956. 
An analysis of project costs(4) indicated that the channel had been stabilized — ‘ 
: ae at a cost of approximately $40 per foot through a reach of channel that varied & 

eS _in gradient from 6 percent in the lower portion to 10 percent inthe upper. 
Thisisa considerable reduction in cost over previous 


eniteria with a few minor modifications. These projects are in Monrovia 
_ Santa Anita Canyons, above the cities of Monrovia and Arcadia, respectively. ae 
ba The Santa Anita Project(5) is of particular interest because it is the first 

- attempt to stabilize all main channel reaches of a somewhat larger watershe: 
miles). This project, when completed, will consist of 70 
_ dams and will stabilize 50,000 feet of channel. . ‘The cost of this project is 
_ estimated to be $1,750, 000. Fig. 4 shows the over-all stabilization plan. 
Fig. 5 shows Dam 4 which is typical of the concrete crib structures of this } te ee 


_ The system will provide over 757, 000 cubic yards of debris storage P 
capacity 1 within the > watershed. This capacity is of special importance b be- 


ready lost 56 percent of its capacity through siltation which has seriously 4 
_ impaired its usefulness for flood-regulation and water conservation. As an 
_ example of the debris movements that take place in the San Gabriel frontal _ 
watersheds when subjected to a fire- flood sequence, two relatively minor 
storms that occurred in in the Santa a Anita were hed & in J anuary 1954, following 
Baas bur rning over between Christmas and New Year’ s, ONS 208, 100 cubic 


i largest of these was the Brown’s Mountain Barrier which rises 58 feet above phe 
a: Be stream bed and has a storage capacity of 1,200,000 cubic yards. A total of a 4 ye 
67 ec were installed using five tyne af conctruction materiale: 
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MOUNTAIN CHANNELS 
y larger of the two storms. During the 30 years ‘since the reservoir was com- 
& pleted, debris has entered at the average rate of 74,200 cubic yards per year. 
4 Assuming this rate to continue, the storage capacity provided by the project — 
is ample to hold all debris eroded within the watershed for more than 8 years. 
<4 Through sluicing finer materials reaching the reservoir into the downstream | 
hae channels, it would then be possible to reduce storage loss to about 35 percent _ 
a. of the original capacity. It is also anticipated that the debris production from _ . 
a the watershed will be greatly reduced and that any debris that does ss p the © a 
reservoir will be fine material which can more readily be sluiced. 
a _ Water conservation benefits, as previously discussed, will be realized trom 34 


the 164 acre-feet of water storage capacity that will be available inthe im- a 


pounded debris. This volu 

a percent. The pore space in the 70 small reservoirs will fill during times of - 
storm runoff and drain out slowly as the storm subsides. _ This cycle may ee 
— itself several times during a single season, depending on storm spac- 
ing and duration. ‘Therefore, a somewhat larger continuous supply will be a. 

available to the spreading grounds in the valley below and waste to the ocean oa : 
times of sizeable flood runoff will be correspondingly : reduced. ‘The in- 
4 


The conservation of water is also the objective of another important chan 
nel project that is his project being put into operation in Monroe 


he 
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-_ == yield, peak runoff, and base flow i study the effect of treatment on a: 
> ’ se flow inasmuch as continuous rainfall and 

| 


California Forest Range Experiment Station of the Se Forest Service, 
is a study to test on whole watersheds the promising an and practical enthois. a 
of increasing water yield by the c control of riparian growth. | The study will 
ales determine how increases in water yield differ with variations in rainfall, 
= soil, geology, topography, and vegetation. ‘The influence of these factors mus 
be determined because only with these basic answers can this type of trea 
ment have widespread application to other watersheds. Other agencies of 
_ government are also cooperating in this project to determine its effects on — 9 
of water, on erosion, and on peak runoff. The Monroe watershed and 
. adjacent Volfe watershed, which is to serve as a control, are located within Te a 
the San Dimas Experimental Forst, and have been calibrated in 
x side slopes. The vegetation removed consisted of an overstory of large trees 5 
_ such as live oak and alder and an understory | of brush. The treatment, when _ a 
= will involve about 80 acres of the 875- -acre watershed. 
* vi It is anticipated that native grasses will occupy the treated area; however, _ 
ae chemical controls of woody regrowth will be required. The clearing <a 
also includes the removal of all floatable debris from the channel areato 
; prevent temporary damming and bulking of storm runoff. It is too early oe 
the study to give valid costs or specific estimates on effects of treatment. co 
is understood that such results will be released by the Experiment Station -_ : 
; ‘The channel treatment projects in the Arroyo Seco, Cooks, Monrovia, and 
pe. Anita Canyons were undertaken jointly by the Los Angeles County Flood 
Control District and the U. S. Forest Service as part of the Los Angeles pay oF 


2 


River Watershed Flood Prevention Program. should also be noted that 


‘Fig. Dem 4, Santa Anita Canyon, near Arcadia, California 
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much | _ been accomplished in regard to fire improvements, vegetative 
a cover improvement, and road ‘slope stabilization as apart of th 


watershed treatment programs. 


channels of a watershed by mechanical means has emerged as the most 

practical m method of ‘accomplishing | debris reduction. Mechanical stabilization 

ae by the use of check dams reduces erosion by halting channel downcutting and 
by ‘reducing | the occurrence ce of bank sloughing. Other benefits also provided 


_ conservation of water through direct storage and increased percolation op-— 
portunity. The Big Santa Anita Project, near Arcadia, which is presently 


under construction, represents the first full- -scale use of this type of treat-— 
ment in this area. The project, when completed, will stabilize 50,000 feet of ae 
to be $1, 000. ; 


Rowe, and R. ‘Nagel is herewith ireatly 


if Debris Reduction Study Group, Debris Potential Los Angeles 

_ County Flood Control District, ~ 

Debris R Reduction | Study Group, ‘Report o on Channel Stabilization 

Angeles County Flood Control District, 1958. 


U. S. Forest Service Photogra h. 
Debris Reduction Study Group, Report on Cooks Canyon Project, Los 


Gq Angeles County Flood Control District, May 1957 — Weare: 

y vd 
Bs, 5. Debris Reduction Study Group, Report on Santa Anita Canyon, Los Angeles" 
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| —s-' This paper has presented the background and some of the details of apro- | © 
gram of research and developments aimed at reducing the volume of debris 
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ELECTRONIC COMPUTERS USED FOR HYDROLOGIC PROBLEMS 
Francis E. Swain,! M. ASCE and Herbert S. Riesbol,? 


The for solving problems conidia: has led Buspen 
cae - Reclamation engineers to the utilization of electronic computers for suc 
studies. Applications: have been made to problems in in water resources 
r in requirements, water utilization, flood hydrology, and sedimentation, Opti- 

of of high accuracy are obtained quickly and economically. 


“ten The need for the rapid solution of a variety of complex hydrologic prob- 4 i 
1 


ems has led Bureau of Reclamation « engineers to the utilization of electronic. cs 
data- ~processing equipment. Hydrologic problems encountered in the planning, 
design, construction and operation of multipurpose water resource pagal + 

are, ini many instances, ideally suited for the use of these machines since Hoes 
Pe such problems often require the solution of complex mathematical i a 
E ships, the reduction of large volumes of data, , the f frequent repetition of a _ ian 
operation, or or the evaluation n of alternative assumptions and criteria. 
a past years, Bureau engineers have obtained sound solutions to many of > 

ot their hydrologic problems by short-cut, , empirical methods. | Continuous _ 


e -drologic research has, however, resulted in the development of many com- 


plex mathematical relationships, and the operation of expanding hydrologic | oh 


4 §2 and data by the usual manual methods because of the greatly increased costs a 
| Bebe. the water and power users and other beneficiaries. Fortunately, the de- aa 
ht  velapliseel during recent years of electronic | data- processing equipment has 
Note: Discussion open until ‘April 1, 1960. ‘To extend the closing date one e month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 2243is | 
fd _ part of the copyrighted Journal of the Hydraulics Division, Proceedings of the f 
American § Society of Civil Vol. 11, No. HY 11, 1959. ‘tins 


i... Head, Electronic Data Processing Section, Div. of atte: 


4 

{ 
te — 
— 
| 
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_ without greatly increasing his requirements for technical ——— 
benefits of utilizing the new equipment ‘may range from the e simple 
. testing of tried and proven empirical relationships, in light of new technical 
ss knowledge and data, to the actual solving of involved mathematical relation- : 
Pe ships. . The conversion of evaporation pan data to obtain gross reservoir a G 
- evaporation is an example of the former benefit, while the mpeynogeer of a 


large volumes | 
‘the latter. "Whenever an empirical relationship is tested and found valid, it” a 


tena’ on the one hand and a full knowledge of the tiaidlagie requirements of 
specific problems on the other hand. With hydrology studies, as with all 4 
problems for which electronic computers are used, it is essential that the 
zs problems be completely and accurately defined before any attempt is made — 
ath. » utilize such equipment in their solutions. Even though a problem is fully 
¥ _ defined, it may be a somewhat difficult procedure to program the desired ae 
‘solution | on an electronic computer. Without such definition, it usually will 
be impossible to program any solution 
9h The engineer who has hydrologic problems to solve must first — his 
od _ problems. However, he also must know the | primary capabilities and limita- — 
_ tions of computing equipment so that he can determine whether or not his _ 
potential applications are indeed proper problems to put on such machines. 


Full integration of problem know-what and machine know-how is essential to . 


q 


the efficient and effective adaptation of studies to machine solution. A proper 
decision to automate can result in optimum s solutions being obtained quickly Pe [a 
and economically. An: improper ‘decision can result in added confusion and 4 
iti is the desire of the Bureau’s engineers trial of 
- electronic data-processing equipment to sample problems covering the entire 
_ scope of hydrologic problems involved in the planning, design, construction, 
bre and operation | of Reclamation pr projects. itt is believed that this procedure will 
a provide a basis of judgment for further and more detailed utilization of the _ q 
= equipment. Such problems may be divided, roughly, into the categories of eee 
, water requirements, water utilization, flood 
and sedimentation. Within these categories, some of the problems which © 
appear to offer the ‘greatest possibilities | for * machine pos wa and which ar 


fact 
| 
dep 
me} 
_ retained and improved as may be indicated. The complex equations and not 
masses of data are utilized as economically acceptable methods of 
— in applicationof jog 
— ~~ . Herein lies the challenge to the engineer in app lesa log 

The application of electronic calculating equipment to hydrologic problems 
general and the introduction of machine computing procedures into anor- a req 
dati 
| 
4 

wat 
4 
4 (the 


| ELECTRONIC COMPUTERS 
ta To extend basic of streamflow, ‘and other similar 
f 


actors and to fill gaps in those records. acs } 
evaluate changes in observed hydrologic data in relation to. 
r physiographic changes, changes in record keepers and equipment, and ex 
_ To ) adjust historic observed hydrologic data toc constant levels of develo i 
To analyze land use and | watershed treatment data in relation to stream- 
flow, precipitation and other hydrologic data so as to determine or. 
not changes in streamflow can be attributed to land treatment. 


fi To forecast seasonal runoff by use of multiple correlations of meteoro- 

To correlate climatic, physiographic, and othér pertinent data with ee 


To analyze data pertaining to domestic, municipal, and industrial water 


requirements in terms of related factors such as location, climate, and 


conduct coordinated operation studies of multiple-purpose water re- bar 
‘source developments that involve various combinations of water resources 

water requirements, and facilities, bea 


tro and data for specific streams or or areas. 


_ To derive flow-duration curves and sediment - rating c' curves from n adjusted — 
and sediment sampling data. 
To estimate sediment-yield rates by correlation of meteorologic, hydro- 
is calculate total ‘sediment load tran: sported by a stream. — 


vas J 


a 


oe Trial solutions os electronic data- processing equipment have been com- “et 

_ pleted on a few of the problems mentioned above and are described in the has 
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ns of for the proposed Waurika and Arbuckle 


in were for estimating water su, supply and 


in ‘neighboring areas of ‘These re 

. lationships were explored by machine methods , by developing a multiple 4 a. 

+e ak. regression of monthly streamflow in relation to concurrent precipitation and 
antecedent moisture. Streamflow records which were utilized are: 
Rock Creek near Dougherty March 1956 to 


ela Beaver Creek near Waurik 
Neighboring areas 


October 1943 to 1957 | 
June 1936 to September 1956 
Little Wichita River near 1941 to December 1955 
eS vy iid The analysis involved a total of 282 trial correlations from which the oe 
consistent and reliable estimating equations were selected. Equations were 
developed first for Little River, Blue River, and Little Wichita River from a 
monthly and seasonal correlations of precipitation and runoff. These ec equa- wee 
; tions were applied to the Rock Creek drainage by substituting Rock Creek poe \ 
oe precipitation data and adjusting for area of drainage basin. The resulting ew 
estimates of streamflow were compared with available records of Rock | = 
ES Creek at Dougherty in selecting the most suitable equation for extending the 
Creek records. Although the streamflow records on Beaver Creek, 
_ Waurika Project, appeared to be adequate for the derivation of a precipitation- 
runoff equation for extending streamflow records, the equations for the other ky 
basins mentioned above were used as a check on the Beaver Creek estima- 
ie Multiple-correlation computations for the 282 equations involved in these 
were made on an electronic computer in about 7 machine hours. At 
«$80 per hour, the machine rental was 9560. _ Labor charges have not yet been mn 
fully defined but will be about $2,500. A direct comparison with a — t355 
_ solution is not available, but previous experience indicates that, for say 


detail, manual solutions would have been much more expensive. Actually, 


3 


"sultant accuracy would have been reduced. 


Use and Watershed te in Relation to Streamflow 


An extensive study of the effect of land use and watershed heanbeecet 00 
o streamflow is being conducted cooperatively by the Agricultural Research ee 
_ Service and Soil Conservation Service of the U. S. Department of Agriculture — 
and the Bureau of Reclamation of the U. S. Department of the Interior. tt is 
hoped that one of the products of this study will be a reliable expression = | 


determining runoff asa function of of | related fi factors. Since there are many 
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a means for obtaining such an expression or estimating equation. 


The study group examined precipitation and runoff records for many 
areas) ranging in size from small experimental plots to large drainage basins. Ae hat 
Exploratory have included various combinations of 


possible. Manual ‘would have been prohibitive due to the expendi- 


ture of both the time and money required. 


Runoff during the snowmelt season, which is of major concern in project — 
planning and operation, is influenced by severa! factors on which data te" 
available. Several of the | principal factors are water content of snow antece-— : 
dent to the runoff season, previous precipitation (as an index of watershed aa 
conditions), and precipitation | during the early part of the runoff period. 
} of these factors involves numerous measurements which require investigation : 


a As described by Ford, (1) “the problem is to pes a multiple- Ss 
. equation which summarizes the relationships of past hydrologic eventsas =~ 
evidenced in records of natural runoff and | of factors | contributing to runoff, 
_ and to determine the . reliability of this equation as a means of forecasting. ” » 
2 In past years the Bureau of Reclamation has developed numerous a 
equations by application of multiple-correlation analysis to many river 
in the western states. Many of these computations have been accomplished 
Any exhaustive study of forecasting possibilities would involve large eathhe: 
amounts of computational work— prohibitive amounts so far as manual 
- cessing is concerned. Extensive machine processing of data for the develop- is a 
_ ment of forecasting procedures for the operation of a Bureau of Reclamation eam 
reservoir system in northern Colorado and southeastern 
_ Wyoming was conducted over a period of 3 years and involved more t an 500 i Aa : 
- trial equations. _ This study has been described by Koelzer and Ford. 2) Con- ~ 
sistent improvements in forecasts were made during the period of the machine _ 


Coordinated studies of systems of multiple- -purpose reservoirs 


ig and powerplants pose one of the ‘more intricate and time- -consuming hydro- 
er logic problems : solved by engineers | s of the Bureau of Reclamation. _ These ory 
ke studies may be prepared to serve a variety of purposes, such as: determina. - 
tion of f reservoir, powerplant, outlet t and conveyance system capacities, ef- a 
‘i fective regulation of streamflow to | provide dependable releases for irrigation, 
5 municipal water, prior rights, hydroelectric power generation and other water 
uses; flood control; or economic analysis. A typical study involves the | utili- 
zation of historical, depleted, or adjusted streamflows at specific points. 
_ within a river basin. _ The flows are routed through a system of reservoirs _ 
powerplants into a conveyance system a and to points of use such as ble blocks 
of irrigated land or municipalities. Losses due to evaporation, seepage, q 


a 
we 
ivi 
basis in an effort to develop rainfall-runoff relationships 
b 4 under changing levels of land use and watershed treatment. 
exhaustive study is continuing. Electronic computers have permitted 
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storage a are taken into account as they occur in the process . The studies are 


Bee prepared to cover an entire period of record of perhaps 25 to 50 years, 


. It may be expensive to plan and | 
program a complex operation study, 2 and a number of studies usually 1 must be — 


— to realize a Saving over manual computation methods. — The studies * 


criteria. The e precise p planning, ‘programming, testing, and coding of 


_ tions can be obtained by changing only certain n independent variables or con- 


the operation studies may be tedious, time consuming, and but 
savings in time and money can be realized for a series of studies. 

‘Recently, the necessity for computing a series of ¢ on | 

_ ‘Upper Missouri River offered an opportunity to utilize electronic data- panies 


Processing equipment. Int this instance, it was desired to operate. the 


fv - powerplants to determine the hydroelectric power production for various con- 
if ditions. More ‘than 4,500 00 individual entries 's were needed to compile the re- _ ell 
quired basic data. . These data included the runoff records of the Missouri 1 
_ River and its tributaries from Townsend to Great Falls, Montana; potential — 
r depletions; evaporation rates; reservoir area-capacity tables; power a 
a a and | pumping plant characteristics, such as efficiency, capability, and capacity; 
runoff forecast equations; and flood-control limitations. About 3,500 locations 
of magnetic drum storage were utilized to store the program + and data, or, in — 
other words, 3,500 separate ‘computer instructions and values of initial input a 
data were required to direct and compute the study. = 


With the development of the Modified Einstein and publica- 
: tion of a step method for computing total sediment load by this procedure, (4) ; 
the Bureau of Reclamation initiated a total- load sampling program on vari- 
ous streams where sedimentation is a critical problem in the planning, design, 
ed _ and operation of Reclamation projects. . The total-load sampling programs on 
the ‘Middle Rio Grande and Lower Colorado River were initiated to aid in — 
eS sar comprehensive channelization plans and designs. On the Lower Colorado oh 
= ‘ River | between Davis and Imperial Dams, 18, 88, 102, and 104 4 total | load = 
a ao ‘made during 1955, 1956, 1957, and 1958, respectively. Since 
ee Modified Einstein Procedure requires long and involved calculations, the as 
decision was made to the possibilities oft making tl the e calculations 
on an electronic computer. In connection with Missouri River problems, the 
Division Engineer of the Corps of Engineers at Omaha, Nebraska, had pro- : 
grammed a total suspended load procedure on an electronic ‘computer. The = 
total suspended load calculation is a a part of the total load calculation by | ;the 
Modified Einstein Procedure, and similar functions are used in both calcula- 
tions. In collaboration with Corps of of Engineers’ personnel, a | computer flow 
chart of the Modified Procedure was developed. The pr program has been com- 
_Pletely coded and checked out with sample problems solved previously by 


I 
when full reservoir capacities willbe utilized. = = = § tudes 
it aa = The electronic data-processing machine is a rapid and efficient device for | migh 
term 
propriate tests being incorporated in the computer program to permit the load 
choice of the proper relationships to fit the given problem conditions. Flexi- annui 
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program was developed primarily for the computation of total lentes in 


the Lower Colorado River. It is, however, a general program which can | _—i 
; used to calculate total loads for any stream. The po possible ranges in magni- ics 
_ tudes of input data and output results, together with foreseeable changes whic 
might be made in the basic computation procedures, were all studied prior to _ 
programming, so that basic data from almost any stream could be used to de- = 
It appears that the development of the computer program for the Modified a. 
| Einstein Procedure will result in substantial benefits. Assuming 100 a 
load sampling measurements per year, the amount of time and money spent _ ass 
~ annually on manual calculations, including checking, would be about 75 man- 
- days and $1,900 , respectively. Initial costs of program development will be are 
- about $2,500 but, after that, annual costs for the computer processing of 100 
- sets of sampling data will be about one computer-day at a cost of approxi- 4 _ 


is not completely satisfactory a at this time, and of 


ry the initial program can be changed easily to investigate modifications and es 


‘om 
«a To athe the basis for a reliable estimate of sediment inflow to Glen 


Canyon Reservoir, which will be formed by Glen eye Dam now being — 


airtel anaaies at poate Canyon and Lees Ferry were used in the study. Th 

Lees Ferry station measures the approximate contribution of the main 

Colorado River to Glen Canyon while the Grand Canyon station <j ae 
its 


contribution of the Paria and Little | Colorado River in addition to the main 


Some preliminary studies were conducted to determine seasonal breaks 
in the sediment rating curve and period breaks in runoff characteristics. _ 

_ From these studies, the yearly runoff variation was divided into two parts | at 
for each station. The he sediment v variation was divided into four parts | for each 


~ 


_ The study was performed at a cost of approximately $1, 500. ‘If the total 


"study had been done by manual computation, the cost would easily have (ex- ~ ag 
ceeded $5,000. ‘One important benefit of this study is that the complete runoff 
_ and sediment records through 1955 are available for any other * type of study ae ers 
_ which would utilize part or all of the station records. With the addition of a : 
_ few years’ records, the study can be kept completely. up to date for a small 
cost. _ Variations of the study are possible, using different periods or season-— 
al breakdowns, with little effort or cost. _ Manual computation of study varia- na 
study was simple in to some on elec- 
tronic ‘computers, ‘but resulted in an analysis of approximately 50 years of 


= for a cost and ina comparatively short 
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cam Approa ach to Computer Solutions tie 
ie When hydrologic problems are solved manually using a desk- ~type calcu~ te 


= — to be made and to summarize the individual numerical operations. As values” 
ee are required for measured data, trigonometric functions, logarithms, powers» 
a of numbers, or quantities determined from particular mathematical relation-_ 
é ships, they are obtained by the engineer by searching tables, reading charts, _ 
. or scanning tabulations of experimental or computed data. The order of mak- 
ing the computations and the procedures used are usually so chosen that the he 
calculations will proceed in a manner as efficient as possible. tere) phil 
‘ -—tydrologi electronic computing equipment is utilized in obtaining solutions for 


- the manual approach. : The order in which the calculations are made is so Fig 
i _— chosen as to take full advantage of the inherent logical capabilities of the bhenit 
_ stored- -program computer; such as the ability (a) to choose the proper com- 
putational steps to be carried out for a particular problem through the use of | 
“branching” operations; (b) to modify instructions according to a prechosen 
plan and the conditions of the problem; and (c) to perform repeated calcula- — 
4 tions of a given type through “looping” operations and subroutines. 
A q ‘It is often more efficient to have the computer calculate the valueofa | 
; om trigonometric function each time it is needed than to store a large number ~maile 


tabular values of the function in the ‘memory of the Such an 


- 
"gree of accuracy by a a mathematical expression obtained to ‘curve- e-fitting” 4 
methods. Here again, the computer would be given instructions to 

Value of the expression each t it is required. When tables of are 

such that they | cannot be expressed properly or or adequately by a 
ts relationship, it is often possible to have the computer determine —o. 
vee values tc to be used in a given calculation | through the use of “table look- ae 

“ commands or routines. When a study requires the u use of a mathematical re- 

_ lationship which is difficult or impossible to express in elementary functions, — 

(5) Bt 
approximate expressions may often be emp oyed with a computer tc to give 
4 the desired results. New techniques and procedures in ‘numerical analysis 
cd are constantly being developed for use with electronic computers which per- 
‘mit the application of these machines in many hitherto con 


Many programming aids are now available which ch help to reduce the amount 


of of time and effort required in t the ie programming and coding of problems for 0 
‘computers. Through the use of assembly and compiler routines, it is posni- 
ble to have the computer itself assist in the preparation of a program for a 
given problem. Such procedures illustrate further the fact that the approach 
a the solution of a problem using an electronic computer may be penta 


i 


‘The use of electronic data-processing equipment on hydrologic problems, i 


such as those noted in this paper, has been studied to determine the as- Fase 
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ELECTRONIC COMPUTER 
- limitations of the machines. , It is apparent at once that, while the Memos 


engineering analysis, experience, « or ‘The ‘most obvious reasons 
_ utilizing such equipment are to reduce the time and cost required for calcula- =“ 
_ tions; to provide more data upon which more intelligent decisions can be ee cm 

ae ‘based; and to free the engineer from the > drudgery | of routine manual computa-— is 
_ tions and permit him to do more important professional work. That savings ~ 
qi in time and money can be substantial has been | demonstrated in previous dis- — 


@ we 


‘From a purely technical viewpoint, one of the major advantages of ma-_ 


_ chine ‘Solution lies in the ability to obtain information beyond the economic Ame 

limitations of manual processes. For a given problem, a variety of criteria 

and basic assumptions can be examined with relatively simple adjustments in 
Use of the electronic computer also p permits a more e detailed solution 


of the planner is the ability to compare of alternative 
schemes of development. Reservoirs, powerplants, conveyance structures _ 
2 and other physical works } may be left in or taken out, and i capacities 'S may be mS 
4 changed with only slight modifications in the computer program ‘once it has oe Uh a 
worked out correctly. The rerun time is usually a modest item. 


f may be made and completely rerun with little increase in cost. The addition at 
< or deletion of major project functions such as irrigation, power, or flood » 
Bier usually requires: a basic change in program and, therefore, may re- 41 4 
sult in substantial increases in cost. Once a complete program has been | pre- Es ‘ 
_ pared, say for a multipurpose hydrologic system, alternate schemes of opera- a 
may | be investigated thoroughly, quickly, and economically. end 
"result is the ability to determine the optimum multipurpose operation ‘solution ra 
rather than merely what might be a satisfactory solution obtained by manual a 

m methods. More intelligent decisions can thus be made because more ee i 
is available to support those decisions. 


og 
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Proceedings of the American Society of Civil Engineers 


DIGITAL COMPUTERS FOR WATER RESOURCES INVESTIGATIONS * © a 


Digital computing equipment is now being used by the U. S. Geological Sur- 


= to analyze published streamflow data and to process data obtained in es : 


in the field on punched paper tape. 4 


| Digital equipment is now being used Geological to 


iis | s analyze published streamflow data and to process data obtained in connection 7 Fd. 
with water - loss studies. Field recording of basic streamflow datain digital 
3 form on paper tape ; is about to begin, and a general-purpose digital computer ie 


2 will be used for the computing and tabulating. Statistical analyses of stream- $$ . 
flow data are now being made that hitherto were considered impractical be- 


aa cause of the tremendous mass of data involved. Water - -loss data are being 7 | 


— 


| processed rapidly and the results are available much sooner than if the data 
handled manually. Trained technical personnel relieved 


Note: _ Discussion open until April 1, 1960. To extend the closing date one ne Coa ia 
za request must be filed with the Executive Secretary, ASCE. - Paper 2244 is an 

‘part of the copyrighted Journal of the Hydraulics Division, Proceedings of the Be a 
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American Society of Civil Engineers, » Vol. 85, No. E HY ¥ 11, November, 1959. 


a. Prepared for presentation at ‘American Society of Civil Engineers ( Con- 


Fort Collins, Colo., , July 1959. 


Engr. U. S. Geological Survey, Denver, Colo. 
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possible use of modern digital computing equipment in hydrologic __ 
ve: 


Meh. _ Many different meteorological and limnological data are obtained i in water-_ 

: Manual data processing is slow and laborious. The computations 
if are of a type that could be made efficiently using a digital computer. _ Studies — rs 
i by” the Geological Survey have been directed principally toward the use of Py 
_modern digital equipment for | processing streamflow records and 


_ ‘The Geological Survey is currently making full -time use of its gene 


“a fourth of the available time ‘paling used for hydrologic studies. Other studies ie 


are being made at Denver, utilizing a similar computer at — Research ie 


a _ Institute ona part- -time basis. 


mt Before any decision could be reached as to whether it would a seminal 


> ei the Geological Survey to obtain a general-purpose digital computer, pre- ‘- 


ze liminary investigations were made to determine the type of studies that might _ 


_ be advantageously made with the computer. The ease with which Gtenstve 

ie. - statistical analyses of streamflow records might be done on such a computer ee 
ws was considered. The Geological Survey has been collecting records of 
_ streamflow for many years, but few analyses of the data had ever I been made. — 

; In order to solve certain problems relating to streamflow, a user of our i sitar 

_ records first had to abstract from published yearly reports hundreds or per- 

_ haps thousands of figures of daily discharge for a given s stream-gaging station 


and then arrange and summarize them in the manner ‘required for his analysis 
2 of the records. Several users of the records might duplicate each other’s 


nalysis almost exactly. The availability of a general-purpose digital com- 


_ puter made possible for the first time the preparation of standard statistical — 


_ analyses on a large scale which could be used as the basis for the majority sabia 


5 a studies involving figures of f daily discharge. Three b basic 


scharges were selected as follow 
~frequency analysis of the daily discharges s asi 


duration each water 1 to Se pt. 30); 


ae - computations are relatively : simple for a digital pee’ 


is, frequency distribution class s limits are c chosen ond the 
“number of daily discharges in each get is tallied. The totals in each clas 


time-consuming task, and ited to mechanization. 
{Statistical analyses are needed to make best use of streamflow data. Many also 
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flow equaled. or exceeded each of the class limits. 
q _ The computation of lowest flow for certain numbers of consecutive days * - 
= is even simpler. A running sum is computed for the number of days | . 
i. under consideration. " Each time a sum is obtained it is compared with the — 
lowest sum for that number of days previously obtained. Ifthe new sumis 
than the previous lowest sum, the new value is substituted. Thenthe 
first value of the series of discharges is subtracted, the next discharge ahead ae g 
of the old series is added, and the comparison with the previous low sum is Bon Zz 
- again made. When the end of the year is reached the lowest sum thus ob- 
"tained is divided by the number of days to obtain the lowest average discharge ‘i 
oe that number of consecutive days. Next a new number of days is substi- oe &g 
- tuted for the old number and the whole process | is repeated. . After | the results Ei 
_ for all the required numbers of consecutive days are obtained, the results for — 
_ that year are read out of of the he computer. High-! flow ‘summaries are computed r 
= It is obvious that the scaaeieantiils is simple, . but the number of data is lenge. 
a The biggest problems encountered were finding an efficient method of tran- _ 
_ scribing the masses of data into a form suitable for direct computer input and 
_ developing procedures for verifying the prepared data. For instance, after a 
4 detailed study of various possible storage media, paper tape was chosen. The | 
a most important reason was that the data consist primarily of one value (dis- 
\ _ charge) for each time period (a day) and the data can always be used in ee 
_ rigorous time sequence. Therefore, mechanical sorting is not needed, and 
e E advantage can be taken of the smaller amount of space required by a unit ess i 
a data on paper tape than on cards. Also the ability of the Datatron 205to — ond - — 
= read paper tape photoelectrically at high ‘Speed (540 characters per ae ; a vs 


a4  -—_aa have been used for most of the paper tape preparation to ae: ) es 7 
_ These machines add and punch simultaneously, which allows using previously * 
computed totals on the source document : as a check. _ However, data for com E 

_ puter input must be practically without error if the computer is to be used 


: rrors in a the. tape, this procedure has proved to be the most practical solu a j 
. tion, and the computer is now used to edit all data tapes. About 25 seconds pst ‘oa 
~iper station- year is needed to to make six validity checks. After the machine 
editing is done, the actual computations, which require about 4- 1/2 to 5 p.gee! Ses 
_ minutes per station-year for the three basic analyses, are seldom interrupted 
because of errors int the e data. To date the 1e computer I has been used to imelpee 
: Much work has also been done in developing means for using the computer Re 4 
2 in processing raw streamflow data. Preliminary studies: were begun long 
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computer application is much more complicated than results 
ready c and a solution is only, just now now at hand. 
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November, 1959 
For y years stream have been recorded on strip charts. This is an 
analog record. As processing these charts is laborious and time- consuming, — 


puter to scan 1 the strip charts | photoelectrically, gage 
height at uniform time intervals, converting each gage height to its corres- __ 
ponding figure of discharge, and integrating to obtain the daily total. ae 
basic design of this special-purpose computer was simple, but many unantici-— 
_ pated complications were encountered when the computer was built. It was 
ey completed, however, | and v was operated experimentally for a short time. Al- te 
A _ though this device could not be made to operate with adequate reliability and | 
_ accuracy for continuous use, a great deal was learned from it. The most im- *9 
portant thing was that special-purpose computers are inherently expensive, 
oa aa because development and testing costs cannot be spread over many similar ee 
: a devices as they can be for general purpose computers built in quantity. 
Therefore it was decided to utilize general-purpose as as 


a If the basic data were recorded directly in suitable seine, a general- purpose 
Be a could then be used for almost all of the data reduction. But until leg. 
_ very recently the available recording instruments were either too expensive | 
4 or required too much electric power. Cost is important because several erat, | 
: _ thousands | of field it installations are involved. Standard electric power is not b's 
aa available at most gaging stations. A slow paper -tape punch especially de-— J 5, 


_ signed to produce punched tape exactly as required by our general-purpose 


tion and cost too much. A new v device developed within the last year seems to a 
<3 be the solution to the problem. This is a simple mechanical coding device at- 
ie. tached to a slow-speed punch which will record a four- -digit observation ina 4 
_ single row of holes on a wide paper tape. As originally designed, the punch ~ Pia 
hada capacity of only three decimal digits, and required 110-volt power, but B | 
_ the manufacturer modified it to record four decimal digits and to use battery 
power. Fortunately, current drain is extremely low, so that a low-cost 
cell battery should operate the device for a year 
this parallel-coded wide paper tape can be used in one commer- ; . 
cial paper-tape reader, it is not suitable for direct input to any of the present — 
general-purpose computers. A special paper-tape translator had to be built | io 
_ to repunch the data contained in the parallel-coded wide tape onto serial-_ role 
; coded tape. Computer instructions are inserted automatically by the transla- 
care _ tor where necessary. However, in one important respect the translator is a - 
bles ssing rather than a burden. Each of the present-day computers using» eal 
aper tape has its own distinctive coding and format for paper tape input. if 3 
thousands of gaging stations were equipped with punches which would produce i 
paper tape for one specific computer, newer computers: ‘could not be used un- 
less each field instrument were modified, which would be prohibitively ex-_ ae 
pensive. By using the parallel-coded wide e tape, the data can be recorded in 
the least possible space and the translator, with its easily changeable coding a 
matrix and sequencer, can be used to put the data into the particular coding oe 
_ and format required for input to any computer. 
‘This new digital water-stage recorder | will involve the following sequence 
of operations. _ Stream gage heights will be punched on wide paper tape at the a 
river bank installation at intervals of 1 hour, 1/2 hour, or 1/4hour. About — 


nce a month the tapes will be removed and sent to a local area office for 
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translation ‘into suitable for direct ‘computer e! translated 
_ tapes will be sent to the computing center and fed to the computer, each with , 
- a manually prepared tape containing the stage- discharge | relation for the gag s 
7 ing station. Daily discharges will be computed and tabulated and at the 2same 
& time stored on magnetic tape for future use. The printed tabulations will be | ie 
a sent back to the field offices which operate the gaging stations so that the ee. 
originally computed results can be examined for possible discrepancies. 
te After any necessary corrections are indicated, the tabulations will be re- | Ree ; 
turned to the computing center so that the information on magnetic tape may 
: be corrected. _ After this, the computer can be used to print tables of daily 
_ discharge in a form suitable for direct offset printing of the annual series of 
Geological Survey water-supply papers. The data on magnetic tape will then 
~ also be available for statistical analyses such as those previously described. y a 
ee ‘This system of data collection and data reduction is now being put into 
operation. The prototype ¢ of the digital water - -stage recorder has been sub- 
- jected to extensive environmental testing and the production models are _ 
ing off the line. . The special paper-tape translator has been tested and found 


ie? satisfactory. The computer programs involved have been written and coi 


“debugged.” The first 100 field installations will be made this summer. _ 
Within a few years we expect to have thousands of such instruments in ie 
field and a full-scale automatic data- -collection and data- reduction i 
Meanwhile, an intermediate system of data reduction has been developed. to 
permit one of our large field offices to use the computer to do some of the gy 
computation of daily discharges that is now done manually ‘Data take 
- manually from the strip chart are read into the computer. | The output is a ee 


Other types of hydrologic ‘computations have been made in Washington 
_ ‘Many more will be made as the capacity of the computation facility is in- | ‘ees 
wane i. Smaller projects already | completed include > computation of of flow in ee 
tidal streams, statistical correlation of physical : factors influencing | runoff, — 
‘recomputation of special discharge measurements to evaluate the effects of 
_ the number 0 of sections taken, and computation of tables of mathematical func- 
4 tions used in ground- -water hydraulics problems. wk ie 


Use ¢ equipment in 


—— the Lake Hetner evaporation studies of 1951 


3 pace from lakes and reservoirs and has investigated methods of determin- 2 
ing evapotranspiration losses from vegetation covered surfaces. In both a me 
_ evaporation and evapotranspiration studies the basic technique has been the — The 
energy-budget method, in which an accounting is kept of all incoming and out- 
going energy, the residual being the energy utilized for 


fect of changes in energy storage is } taken into account. 


3 
— 
= 
a 
se 
—| 
i 
tO 
q 
im 
Ds 
TLL Ul Tal ine ite lor the solution of the 
~=energy-budget equation can be obtained at a centrally located meteorological 
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the reflected solar 


the incoming long- wave ion 
e reflected long- wave seh ox 


Bien the long-wave radiation emitted d by | the ‘body of water or the 
Ta , the air air ‘temperature 
wet-bulb temperature j= ay 
_ For reservoir evaporation studies it is simpler and sufficiently accurate . € 


Pe compute Q using an empirical method developed by Anderson. (1) For evapo- 


transpiration studies, ‘however, little is known of the ‘reflectivity of the sur - ibs 


face, and it must be measured. For evaporation studies it simpler to 
compute Qbs, as it depends only on the | emissivity of of f water, which has been 
accurately determined, and the water surface temperature, which can be | 
measured with little difficulty. A total radiometer can be 
to measure the incoming radiation. For evapotranspiration studies, 
es however, both the emissivity of the vegetation covered surface and its | tem- + 
ne ies perature are unknown, and it is therefore necessary to use a net exchange 
as radiometer, which measures the difference between incoming | and outgoing a 


- padiation. The air and wet-bulb temperatures are used in computing the _ wr 


eek, _ The output of € each of these instruments is a potential difference or voltage. 
- With the use of suitable » scaling devices , all items can be recorded on one 

multi- channel recording potentiometer. Ordinarily an 8-channel instrument 
is adequate. With t this s instrument, the > Various items may be recorded : at ee 
ey intervals so that with an 8-channel recorder, each item is recorded 
_ Different symbols, such as circles and crosses, and dif- = 
- ferent colored inks are used to identify the items recorded. bs 
For computational purposes, hourly average values of each item are | 
usually required. With some items, such as air temperature or wet- bulb 
temperature, the diurnal change is reasonably regular, and hourly averages ki 
be determined quickly and accurately using a graphical technique. For 
_ the computation of solar and atmospheric radiation, three items must be re- 
nee corded, namely: (1) a voltage proportional to the solar radiation received by 
the pyrheliometer, (2) a voltage » proportional to the total radiation cof the 
by the radiometer, and (3) a voltage proportional to the temperature of eal d bs 
radiometer plate. On clear days and on some totally overcast days — 


~ diurnal variation in these | three items is also regular, and hourly averages” ve. 


_ Can be determined quickly. On partly cloudy days, however, transient cloud — 
“effects cause rapid, short period variations of considerable magnitude 
each item. If hourly averages are | to be determined graphically, it is usually 


4 necessary to draw lines connecting consecutive points. Even so, the con- 0 , 
st structed trace is so erratic that the accuracy of the graphical average ~ ae 


‘The frequency of partly cloudy days is sufficient to make manual data 
processing tedious and time-consuming. Experience has shown that in one 


im 


Be: ‘of consecutive points in a desk calculator and compute the arithmetic average. 
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eg field. The fact that the work is tedious should not be be disregarded sii Piaget) 
: 3 tasks of this nature should be kept at a minimum to obtain maximum produc- 
q _ tion from competent technical personnel. As the operations are repetitive “?. ) 
_ and the field data can be put easily in the form of a digital record, the use of Be 
- modern high- -speed computing equipment was clearly indicated. Because we 
computer was available in the Washington office of the Geological Survey, it Q a. 
was decided to proceed | with development of equipment that would record the rs 
data on punched paper tape | in the field. The possibility o of punching the tapes — 7 
_ in the office was considered, but discarded, for the time required to oie 
- _ the data from t the recorder charts to tape would be as great as hitherto ed 
needed to process the | charts manually. It was considered inadvisable to meet 
eliminate the chart record. Visual inspection of that record often 
_ short periods of questionable data resulting from equipment malfunctions. — cee a 
__ Inspection of the punched tape would not suffice, | and it was therefore deemed 
essential that the data be recorded both on the chart and paper tape. iii waite: 
- i field punching equipment consisted of three components, namely: (1) = 
_ analog- -to- digital converter or encoder to convert the angular rotation of nalal 


(2) a programmer to insert certain computer commands and | 
| lated information, and (3) a motorized punch to record the data and eg i 


first component, the encoder, was commercially available. The device 
divides the 11-inch v width of the recorder chart into 889 finite parts and con- 
_ verts the information into binary-decimal code. 
i eer The second component, the programmer, was custom built. It would tave 
_ been possible to punch only the basic data, a three-digit t number giving an out- 
put voltage. The tape would then have to be repunched completely in the office be 
Ls in order to put the tape in a form usable in the computer. For example, aa 
_ the Datatron 205 computer the data would have to be organized into (10-digit — q 
“words.” Each channel requires identification. Programming instructions 
_ must be provided to identify the end of each hour and the end of each calendar — a) 
day. . The computational program for the meteorological data previously a, 
- named is rather complicated; it was originally estimated that approximately — 


3 two-thirds of the total storage capacity of the Datatron 205 would have tobe = 


wd reserved for the program commands, leaving only one-third of the storage __ 

it i for data. This would have been insufficient for a full day’s data, so that the % 
a , computer was ordered to read in data for a 12-hour period, make the compu- ae 
a i tations, and deliver the results before starting another 12-hour period. Later Re 
_ § it was found that the program required much less storage space than had Pay ea 
been anticipated, and the schedule has been modified accordingly. 


a The programmer is necessarily custom built if the punched tapes are to Roc 
z be used directly in the computer without resort to extensive tape preparation. be 
_ However, it does not seem wise to construct a large | number of programmers 
that will produce tapes usable in only one computer. Computer designis not 
static; tremendous advances have been made in recent years and more are to = 
be expected. The expense of repunching the tape in the office prior to com- | . 


= 


ik puting 3 must be weighed against the cost of modifying a number of field instal- ay 


= the third component, a commercially available punch is used. The only ip &§ 
| required t to make the punch suitable for unattended operation 


> 

— 

| 

a 

¥ 
— 

| 

ki 

— 

a 


an evapotranspiration research project in Nebraska in 1958 field 


were recorded on three multi- recording ah and at 


found that it took approximately 1 hour to compute one week’ s data, as a | 
to 2-1/2 work days required. for manual data processing. 
_ Some possible modifications of the computing program are under vitae ee. 
_ sideration. The actual time required for the computer to make certain com- 
- Putations seems to be rather long. _ Preliminary studies indicate that it might | 
be advisable to compute the radiation items from daily averages rather than 
from hourly averages. Over the diurnal range usually encountered, many of 
the relationships are nearly linear, and the error resulting from using daily — 
_ averages appears to be negligible. A critical review of the complete program jf 
"indicates that certain revisions will result ina substantial reduction in com- -— : 
"puting time. "Modifications in the recording and computing scheme have al- 


The advantages of using digital computers for hydrologic studies are many. 


eg oe In studies using streamflow records, the e computations are not particularly 
; complicated, but the mass of data is tremendous. - The use of digital comput- 
‘a io ing equipment will permit much better utilization for statistical studies of th : 


a ne not as) many y data are involved, ‘but the computations ' may be Selinten 
ee - more complicated. For all types of hydrologic data obtained by the Geologi- 
$c, cal Survey, the trend is toward obtaining a a record in the field on tape which — 
can be used directly or automatically converted for use in a digital pc pl 2 
: itt ~ With high-speed computers the results can be available much sooner than 
f 


the computations are made manually. For Various | reasons 
4 


_ With digital computing equipment it might be practical to increase the num . 
ber of such records substantially, for with this equipment the cost of recom- - 
- putation at the end of the water year (which is almost always necessary) — 
Another advantage resulting from the use of digital computing equipment — : 
or hydrologic research is the ease with which the basic | data c can be placed — 
a in a form suitable for publication. The basic data for some research projects 
ee have not been published simply because of the labor involved in preparing the - 
tables and setting them in n type. . Modern computing | equipment can produce pa 
tabulations of data and results of computations that are suitable for offset oie 
reproduction. Column headings and rulings can be added easily using me 368 
Perhaps the principal advantage of modern computing equipment faccrues 
from more efficient utilization of technical personnel. Much of the data 
processing and many of the ‘computations required for hydrologic studies are’ 
both monotonous and tedious. — - Reducing the amount of time that trained per- | 
; sonnel must spend on purely routine tasks will pay dividends in more com- 
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HYDRAULIC DOWNPULL FORCES ‘ON I ‘HIGH HEAD GATES# 


gp 


‘The pay paper discusses model studies in the Bureau of Reclamation’ Hydrau- 
lic Laboratory in which both the direct weighing method and the pressure area 
computation n ‘method were used to determine hydraulic: downpull on two > high : 
head gates, and one model “study in which the pressure area computation 
method alone was employed. . Also described are the methods used and 
4 


_ The general design of a gate hoist provides for sufficient capacity to open 


_ or close the gate under full operating head. The required capacity of the es: M 


. hoist is usually obtained by totaling the frictional resistance of the gate, the __ ee i; 


é “dead” weight of the gate parts being moved, and the hydraulic | forces }tending 4 
ae bas to open or close the gate. The weight of the gate and other moving parts, eae id 
“§ either dry or submerged, is easily computed. The frictional resistance can 

computed by using assumed coefficients of friction. However, the hydrau- 

by Ghe bove and below the gate 
o§ ic force caused by the unbalanced water pressure above and below the ga a 
deaf particularly when the gate is in the throttling position, is far from re 


_ simple to determine. A few tests to deter mine these hydraulic forces are ‘ 


passing beneath a gate leaf creates a pressure due to a 

* 3 _ Change i in the direction of flow and a change in energy from pressure head to 4 


velocity head. This reduced pressure ; acts on the bottom of the gate. There 
i is a pressure on the top of the gate, water pressure if the leaf is ie teem ps 


4 


| 

— 

Note: Discussion open until April 1, 1960, To extend the closing date one month, a 
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fei hydraulic | force on the gate is equal to the difference in pressure above an 
4 _ below the gate in pounds per s uare inch multiplied by the affected cross ae 
gat qua 


i ‘sectional area of the gate in square inches. By changing the shape | of th 
om “ _ bottom of the gate, or by altering the flow passage, a pressure change might 


at > installation offers a new and specific arr ag in regard to hydraulic sy 


Fixed- -wheel 


studies have been mate to determine the downpull on 
several different types of gate installations. One: of the Bureau of Reclama- 


- tion’s earliest model studies concerned an . emergency bulkhead- -type gate for 
Shasta Dam, designed to close the inlet to a conduit discharging under 323 ba 


was 


_ suspending the gate froma spring scales and taking continuous weight read- othe 

_ ings while the gate was being raised or lowered (Fig. 1). However, the hy- oe 

-draulic downpull values were extremely difficult to determine because vol 
excessive and variable frictional forces. The roller chains were removed ne 

_ from the initial model gate a and replaced with small flanged wheels resembling 


i head . An attempt w was made to “weigh” the model downpull forces by cr 


_ miniature railroad wheels. _ Under heads up to 10 feet the apparatus omen 


, too operate satisfactorily, but for higher heads the wheels would stick momen 
tarily and release suddenly, producing the same untenable conditions en- 
countered with th the roller ‘Chains. A new model gate was fabricated using — 


- wheels having a a diameter eq equal to the thickness of the vr: oe wheels a 


Piano wit 


a achieved which would reduce the hydraulic downpull. However, each new anal 
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Vs tremely difficult to record. It was found that the roller tracks had to be imide se 
glassy smooth before acceptable readings could be made. Sey 
4 a Several strategically placed piezometer taps were installed in the model — Pe 
gate and the downpull forces were computed by the pressure differential-area 

method. Results of the for various gate 
_ From these tests, it sgpented’ that the extreme precision required in 

- fabricating this type of a model gate to be acceptable for weighing the at. Po 

forces overshadowed the tedious recording and computations required in 
_ The gate developed in the above studies was installed at Shasta Dam, ad 

- an opportunity arose to measure the downpull forces in the field installation. _ oe 


a _ An SR-4-type strain gage was bonded to the gate hoisting stem and the change - 
= 4 in resistance of the strain gage was measured by a portable strain indicator. 7 


gate position was measured with an engineer’ chain. As the test 
es ceeded, the strain indicator, the engineer’s chain, and a timing device were» a: 
photographed by a motion | picture camera operating 3 frames per second. 
4 _ The timing ‘device was included to correlate the | gate opening and strain read- be 
_ ings with a metering device. _ A cchart showing the hydraulic downpull deter - 
mined by a model study using the pressure area method a and by field meas - buns : - 
urements is hown in Fig. 3. The discrepancy between ‘model and prototype 


at about 20 percent gate opening is perhaps due to improperly computing the | 


. point NY 300 WHA 
+5200 
FORCE = 150,000.89 “rons i000 8s 
i = | FORCE: 


BOTTOM GATE WITH EXTENDED LIP 
| 
SHASTA DAM OUTLET COASTER GA 


PRESS! 


PRESSURE STUDIES OF GATE BOTTOM 


SLOPING | BOTTOM GATE 


— 
yn 
turned on ball bearings to minimize the rolling and bearing friction. Never- $7 ha 3 
: 
if 
x, 
ra GN 5 
— - 


eben, 4 1959 


downpull contributed upper seal ‘it er 

‘Hoo 


_ A laboratory study was made of | the ubdradtic downpull forces on the — 
‘ 4 cylinder- type gate at Hoover Dam. Since the horizontal forces acted a 


— and were resisted internally by the gate, bearing and rolling friction 


the simplest method of solution. 
‘The cylinder gates were ‘suspended by molybdenum - wires in a 1:24 
4; seca model of one of the intake towers (Fig. 4). The change in downward 


forces ¢ on the | gate oe be er from the measured change in strai 

_ At openings of 20 inches and larger, t the reduced pressures on the bottom — 


w tic 
range, a small change in gate opening was “accompanied by only a a slight 
= change in the hydraulic downpull, and the model apparatus was quite stable, vie 
producing reliable results. AS the gate 1 neared the closed position, area’ 
of the flow p passage between the surface on the bottom of the gate and the gate 


a i seat increased in the direction of flow producing a venturi effect which tended | 


tof further reduce the pressure on the gate bottom. $ At these smaller gate oa 


it the gate between the guides would change the shape of the flow passage sander 


La These } pressure changes would be either uniform around the periphery of the 
Ae gate or unbalanced, and the gate would tend to move in the direction of the © 
lowest pressure. The gate movement wo would I resisted by the 


of | the gate were due + mainly to the change in direction of the flow. - In ae 4 
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ice: a tribution of forces would shift rapidly to a new location, and the entire model, § 2a 
although sturdily constructed, would shake so violently that the gages could 
“not be read. No practical solution to the problem was found, although many 
lm MEN: 
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ve corrective measures were tried. For t the final tests, the poe was removed, 
; piezometers installed, and pressure area studies were made with the e gate 
wedged tightly in predetermined positions. 
: is In the field, downpull tests were made on the Hoover Dam cylinder gate for 
ns openings up to 15 inches. Prototype downpull was measured by means ae 
ee strain gages ‘on each of the three lifting nuts (Fig. 5). . Continuous oscillo- | “Ee 
a graph records of the change in resistance of each gage wees made during the " 
tests, and from these records the downpull was determined. The downpull © S 
curve determined by these measurements followed the trend of the model 


4 termined by model study, and the peak downpull occurred at 4-inch 1 gate Bain 
ing whereas the ‘model indicated a peak at 6-inch gate opening (Fig. 6). ee 
Obviously, the of downpull forces on prototype gates 


covering entrance, 20.0 feet. 39.3 feet high. 
closed gate at this location subjected to a head of 156 feet, would be required 
. to withstand a hydrostatic force of 7,980 kips. In another plan, a beam aa 


_ placed across the opening above the ' tunnel entrance to form the top part of a 


_ the semibellmouth (Fig. 7). This arrangement permitted the gate to be in- 
: stalled a short distance downstream from the tunnel entrance where the open- 
ing was 19.7 feet wide by 28.0 feet high. | The hydrostatic force on the closed : 
gate at this locatian would be 5,540 kips, or a load reduction of about 30 per- 
cent. The latter design was preferable economically, and a model study was 
‘instigated to determine the hydraulic characteristics of the installation, par- 


hod 


of determining mga downpull was not favorably considered for the — 


thes obtained, 36 piezometer taps were placed in the 1: 39 scale model gate 
_ (Fig. 8). The approach to the gate, the beam forming the top part of — 
semibellmouth entrance, and the transition downstream were formed of light- 4 
weight sheet metal (Fig. 9); 24 piezometer taps we were placed in critical points — 
the beam, entrance, ‘gate slots, and transition. entire apparatus was 


were e carried habia ' the wall of the tank and attached to a gage board. i 
Preliminary testing concerned the general operation of the system to saga 
tain the overall design would be hydraulically acceptable. ‘When the ; ae 


was remedied by ‘making a fairly large recess in 


- downstream face of the concrete beam (Fig. 7), thereby keeping a balanced — 


| 
— 
a, ___ urements must be made with such extremely confining tolerances that the ey | 
One of the preliminary plans for both the power tunnel and the outlet works 
— 
— 
— 
— 
— 
| 
— pressure areas might occur on the gate so that an accurate pressure study = aaa cs 
= 
— 
pressure on the downstrea ce te exceeded that on the upstre - 
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oe discharge through | the penstock was s computed to to be 4 45, 000 cfs. ae 
_ Preliminary downpull t tests disclosed that the forces prevailing | with this « dis-— 
; charge, and with the gate fully opened, were about as expected, producing hy 

draulic downpull forces of a little more than 700,000 pounds - However, a — 


why A study o! of the pressure distribution on the gate bottom disclosed that con 
_ siderable variation in pressures occurred between adjacent piezometer taps - 
In computing the downpull, therefore, the small areas adjacent to each of the — 
45 taps were considered separately. The sum of these individual downwar: | 
forces, together with the downpull attributed to the upper and lower seals, _— 


constituted the t total hydraulic downpull for any given condition gate Opening 


Model results 


200 300 400 
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rk ‘open, and the gate then ained the hydraulic downpull would vary as shown vin 

_ ‘Fig. 10. . The maximum hydraulic downpull of of 895,000 pounds at about 5 

4 percent gate opening was peculiar to this particular g: gate bottom sinticicehion ba 
however, the 707,000-pound hydraulic downpull at 100 percent gate opening 

7 would exist for any gate with the same cross 5-sectional area regardless of gate 


a An extension of the downstream lip of the gate would create a ogni agi 


yet y 


"GATE SECTION, 


FIXED WHEEL GATE 


GENERAL INSTALLATION 


1 


November, 1950 


é 


‘The maximum hydraulic downpull for a 29.5-inch lip extension was 895, 000 


: i pounds, while that for a 68-1/2-inch extension was 670,000 pounds. ae | 
added cost of strengthening the 68-1/2-inch lip to withstand a thrust of tis) Tok 


; he a 270,000 pounds (Fig. 10) was greater than the ‘Savings which would be phil 


es realized \ using a smaller capacity hoisting cylinder. r. The designers, there. ee 


_ fore, placed a limit of 550,000 pounds thrust on the extension, or a lip ex- bio: 


ia ‘Some reduction in downpull could be a achieved by reducing the length of, uaa 
an 


d streamlining the upstream lip of the gate. The amount of streamlining : 
shortening was again a ‘Structural problem, so computations were made and ‘ 


za the uy upstream lip changed accordingly. The configuration shown in Fig. 11, 
Design No. 2, was considered structurally sound and is the change tested | in 


as _ The fixed-wheel gate would be required to stop the flow caused by some — 


: - operative in the opened position, etc; 45,000 cfs is the maximum possible yet. 
eae flow, and an air vent downstream from ‘the gate limits the head here to 15 


feet below atmospheric. The maximum reservoir elevation is 178.9 feet 
above the lower gate seat. Realizing these e limitations, a calibration chart 
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te openings. Then, for several Res. 


discharges, and each 5 percent increment of opening. ‘Next, by 
_ plotting, a chart was drawn showing the hydraulic downpull which would be | ae 
encountered as the gate was lowered from fully opened to closed, , for any ini- 7 
“tial « emergency discharge. AS the gate was lowered into the stream, the ais 
~ charge and the head at th the air vent would uld decrease until the i the head at the sie : a 
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wee 15 feet. At this ‘maximum m downpull 
: In the downpull chart, Fig. 12, the gate opening is based on a total trave 4 
of 29 feet; however, the gate does not affect the flow during the last foot or so 


The maximum hydraulic is 813,000 pounds and occurs with 
an as the gate 57 percent open. The reduction in the upstream lip extension of 3 
a 11/16 inches, comparable to extending the downstream lip this =" 4 
#3 - would have reduced the maximum downpull by 40,000 pounds. . The 9-inch ‘ay 
ie radius on the upstream lip evidently reduced the maximum downpull by 42,000 
As a result of the hydraulic downpull tests, the hoisting cylinder and con-— 3 
. necters were made sufficiently rugged to lower (or hoist) a total of 1, 300, 000 oh 
A comprehensive field test concerning the hydraulic downpull at Palisades 
_ Dam has not been made. However, during a field operation for another pur 

; a pose, the hoisting cylinder oil pressure was measured while the gate was 
moved from 100 percent open to 89.4 percent, and back to full open. The flow | 

_ through the penstock at this time was 6,000 cfs. . The results computed from _ 

this very limited study are shown in Fig. 12. ‘The trend follows the model © on 
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GATE SEALS AND GATE BOTTOM SHAPES 


have been defined to a degree that will permit an accurate 
ie analysis of each individual problem. Scale model investigations are of pri- a 


: _ mary importance to the engineer in the solution of these problems. The two 
common procedures employed in a model are the direct force 
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must be fabricated to extremely close tolerances. This type of model can be _ 

== investigations should be made whenever possible, and the results pe 
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free discharge conditions the hollow- ~jet valve has proved to bea 
satisfactory control device for flows through large cc conduits operating under 
_ high heads. Results of field tests with a 96-inch hollow-jet valve have re- te Na 
vealed close agreement with hydraulic characteristics predicted from model 
studies. Piezometric measurements, thrust determinations on the valve 
and rates of discharge were included in both field and laboratory 
“as used. The only cavitation erosion ‘evident in the prototype valve was 4 
= byl local irregularities in the body casting, wh which have been eee in sub- 


‘sequent valves by careful foundry practice and inspection. 


flow facta large conduits discharging under high heads, has been in progress 
for nearly half a century. The need for such a control device was responsi- ie 
for r development « of the | Ensign valve (Arrowrock Dam), the needle valve 
- (Alcova Dam), and the tube valve (lower outlets through Shasta Dam). Other 

Valves have also been developed for the same purpose, but the ones named 

a constitute those primarily | used by | the Bureau of Reclamation. All of these 

valves had certain performance and economic limitations. § 

The increasing demand for closer control of releases } through modern 

fs ‘multiple-purpose structures prompted a continuation of studies to obtain a 
mom suitable valve and led to development of the hollow-jet valve. This foots 


type, however, is limited to use as a =“ biti. e valve reventi its 
e one month, a 


: 
DRAULIC CHARACTERISTICS OF HOLLOW-JET VALVES 
] 
va 
«= 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of 
_§@ American Society of Civil Engineers, Vol. 85, No. HY 11, November, 1959. a ae 
_ 1. Engr., Dept. of Public Works, U. S. Navy, Great Lakes, 
2. Hydr. Engr., Div. of Eng. Labs., Bureau of Reclamation, Denver, Colo. 


san of a 12-inch-diameter air model, a 6-inch-diameter xieiealie model 
_ tested in the Hydraulic Laboratory of the Bureau of Reclamation at Denver, es a 
Ta and a 24-inch-diameter model tested at Hoover Dam under a high head. ‘hie. 
The prime purpose of the 24-inch model was to ascertain the hydraulic 
= _ Characteristics of a hollow-jet valve constructed under prototype conditions. 


‘That is, ‘the outer ‘shell, the supporting vanes, and the cylinder containing the 
he. _ needle were cast in one piece. Machine-finished surfaces were limited to the 


ig a ‘Semmes the rough finish of the | casting could ‘conceivably ¢ affect boundary flow 
ys e sufficiently to cause local areas of low pressures, all surfaces in the 6-inch 
The secondary purpose for studies with the 24-inch model was exploration — 
of some critical areas that were too small in the 6- inch model for explora- 


ae by piezometer orifices. A few revisions were ere found | necessary a as a are- — 


Although hollow - jet valves installed at a number of structures, 
initial installation of large units was on the four river outlets through Friant | 
- Dam, Fig. 1. Details of this installation are shown on Fig. 2. One of the four 

96- -inch valves 1 was equipped with piezometer orifices to permit performance 


of special tests to ascertain if this type ' valve possessed predicted hydraulic © 


3 Where dis distances between orifices were small in the direction of flow, the 


% orifices were offset laterally. A photograph taken looking upstream : at this 
valve is shown on Fig. 4. _ Performance of the large valve could conceivably a. 
— from that predicted by the | 6- inch and 24- inch | hydraulic models ¢ due ae, 24 


; A field | testing gg igo sabe: inaugurated about a year after the valves were 

oe valve openings of 10, 20, 40, 60, 80, and 100 percent of needle travel with © 
operating heads of 102, . 180, and 223 feet. . The maximum design head at this 
se installation is 246 feet. For each test, discharge through the valve was ob- : 
a tained from operating records based on current meter measurements in the 
river channel a short distance downstream from ‘Friant Dam. In this instance, 
all flow in the river passed through the valve under test. Hence, measure- _ 
ments of river flow were not subject to possible errors due to subtracting ae 
flow from another source such as a powerhouse 


; Other field observations ‘included general c characteristics of the hollow- -jet 
alves; such as, noise intensity and vibration, inspection of the inte ior of aa 


»sed conduits. The hollow-jet valve ights re- © 
application in closed conduits. tented (Patent No. 2,297,082) with rig 
Bureau of Reclamation and is pa rnment without payment of royalties. = 
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Results of Pressure Measu: 

Pressure measurements were obtained piezometer orifices 


OA valve on each neni line per- 


&g 


-Head piezometer_ 


/<—Direction of 
travel to close valve 


rs yNeedle 


a4 re 3 - Piezometer locations in hol val 


‘Friant Dam Outlet 


pe 4 -Two-rack conduit (horizontal) encasing piezometer 
BE hae leads from interior of 96-inch hollow-jet valve-- of 
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deviatio 
a inlet diameter upstream from the valve and referred to | the same eleva- na fices, il 
_ tion as the other pressures. Hence, the prototype valve was consideredto § ference 
have a horizontal center line that corresponded to the hydraulic models, 
while actually the center line of the prototype valve sloped as | fre AS p 
‘Pressures were plotted against percent valve opening by utilizing a pres- _reveale 
sure factor, F, defined as the ratio of the measured piezometric head, in feet § the valv 
areas 
lished 


ee ——- in the valve by selecting from plotted pressure curves the correct 
Value of and multiplying it the tota total | design on on the valve one diame- 


an example, to find the oh the piezometers, for a 
plot is shown on Fig. 5, when the total design head is 200 feet of water and “4 
the valve | is 50 percent open, , follow the 50 percent line on the plot until it a 
intersects the curve for the particular piezometer and read the value of the | 
é pressure factor at the left. Then multiply the pressure factor by | 200 to ‘ia 
tain the pressure at the piezometer for the case being considered. . If the — 
piezometric pressure is below atmospheric, then the F value is negative ar 
the calculated pressure is also negative. 
Plots o of pressure factors versus valve openings, similar to Fig. 5, were q 
: de for all piezometer orifices installed in the 96-inch valve and in the a 
models. The locations of Piezometers P6 and P25, for which pressure plot 
are shown on Fig. 5, are shown on Fig. 3. 
a ‘The pressure at the needle piezometer in the air space just u upstream from 
‘the vanes was found to be a negative 1.22 feet of water in the 24-inch model 
a when 100 percent open u under a total head of 196.6 feet. In the prototype, the 
corresponding pressure when referred to the same total head was a negative 
4.9 feet of water based on pressure factors obtained from tests at the two & 


highest: eservoir heads. B Based on the pressure factor from test at t the 


variation in the itself, is to the fact that this region ‘ 
is filled with an air-water mixture due to insufflation of the jet in the pecto- 


an increase pressure, while in the ‘model valve, 

little, if any, insufflation occurred. 
Bi. "The ¢ only subatmospheric pressures predicted from model studies were on 
_ the large 1 vanes, but these were not considered sufficient to produce cavitation 
erosion. ‘This conclusion was verified by the prototype tests which revealed § 
an average pressure of minus 18 feet of water in this region, but no pressures 

conducive to cavitation occurred. _ However, negative pressures were found in ‘- 


‘models The average deviation between model and were 
- found to be less than 10 feet of water at maximum Prototyp test head. 


“4 

— 

— 

— 

— : with _ 
total head. All of these values were obtaine 
referred to the same 
— 
— 

— itudes of such pressures were smail and i1 h differed 
res measured in the 96-inc pr 
Ingeneral, pressures 34 h model by an amount approximately equal ot 

those measured in the 2 


HOLLOW- JET VALVES 


‘cause cavitation erosion. of the tested prototype waite 


revealed that cavitation erosion had occurred o on n several localized areas of me. 
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| 


areas was not severe, the metal had been pitted. . Since there was no estab- 
lished or zone of damage, ‘this dition was attributed 


fices, interference by a DO bation of piezometer ori-_| 
terence ead near a plezometer Orifice, a slight dif- 
7 
| 


‘The 
not det 
“held within “specified limits. on the | 
At the time of inspection, the valve had operated for a total of 5,896 hours ‘a the 
at openings varying from 2 to 66 percent under heads from 78 to. 207 fot. BR. 


a Most of the operating time had been at openings less than 30 percent and = ee Pr 


Results of Thrust Measurements 


pressure into the interior, thereby balancing pressures so as to minimize 
power required t to open and close the valve. - Fig. 6 shows thrust on the needle a : 
in the upstream : and downstream directions predicted from the 24-inch model, — ~~ 
te comparable information obtained from the 96-inch prototype. 
For comparison purposes, the units on Fig. 6 have reduced to those 
: applicable to a 1-foot- -diameter valve under a 1-foot head. This permits = 4 
; Bs ce of thrust forces for other size hollow-jet valves operating under 
_ Various heads. ‘These computations can be made by multiplying results” shown 
on Fig. 6 by both the desired head and the square of the diameter of the valve 
under consideration, these values being in feet and square feet, respectively. 
The results of similar data obtained on the 6-inch r model are not shown : since 
locations of the balancing ports established by tests on this small model were 
changed after analysis of results from the 24-inch valve. = = = 
ae ‘The values ‘shown on the plot reveal very close agreement between thrust | 
forces predicted from the 24-inch model and those determined by field meas- 
urements on the 96-inch valve. _ The greatest difference between model and 4 
ee. occurs in | downstream | thrust at a valve opening of 10 per-— 
cent where the prototype value is approximately 94 percent of that determined 
in the model study. The unbalanced force or the difference between the up- — 
stream and the downstream thrusts is the most important, and the maximum 
. ~ unbalance occurs at a valve opening of 100 percent. > This unbalanced force is 
1.29 times the value predicted from the model. = | 


_ ‘The thrust on the needle in the downstream direction was computed from — cs a 


pressure 
= length alon 


radius to piezometer iat 


— 
— 
— 
— 
— 
= 
| 
i 
Total thrust in x-direction = 2, 


oth 
determinable analytically. Values of determined with 
“on the needle were plotted against r values for a valve opening of 10 xra ad 
and the eres. under the curve was attained with a to obtain 

Prdr that value of area ig multiplied by to obtain total thrust 


“the needle in the downstream direction. “The same eget nes was utilized for 
yalve openings of 20, 40, 60, 80, and 100 percent. 208 
‘Thrust in the upstream direction is the pressure inside the needle 
multiplied by the area over which the pressure acts. 
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can result in saving large amounts of money and time. The: time and expense 
involved in performing a field calibration are demonstrated by the fact that Ag 


~ determine discharges through the hollow-jet valves at the hendworke of spd 


_ Friant- ‘Kern Canal, and also through needle valves at the headworks « of Madera 
a Canal. Results show that the current meter measurements could have been — 


aug 
accurate as curves determined by field measurements. 


Once a laboratory calibration has been made of a valve, same calibra- 


ge 


ere e complicated approach ‘conditions disrupt flow 
‘Figs. 7 and 8 present data to support accuracy of the laboratory calibra- — : 
tion curves. For valve openings greater than 15 percent, the variation be- ; 


ae tween model and prototype discharges may be considered as 3 percent. -For~ 
smaller valve openings the difference is greater and may be partly accounted 


to accurate measurements with current meters. 
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MPARISON OF DISCHARGE BETWEEN AND PROTOTY! 


§ oon Comparison of rates of discharge of a prototype and model valve are par- hav 
& ticularly important since they serve to evaluate discharge curves 
— 
: 7 4 Nee years to determine discharges through the river outlet valves at Friant Dam. Thi 
sip 
ytoin- | al 
4 Ulat PUSIUON accurately Peveal true Valve Openings. The 
Particular valves described in this paper are equipped with verniers on the 
— ae 
— 


JET = 


to accurate setting of v valve An equally 
_ important item is that the true head on on the valve r must be known. Difficulties 
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Behavior of the Valve 
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DISCHARGE FORMULA FOR STRAIGHT ALLUVIAL 


A Hein- Kuan Liu, ly M. ASCE and Shoi-Yean Hwang, - A. eli 


ABSTRACT: ‘ati 
a "The proposed discharge formula « contains a discharge coefficient and ex- Bah 
, ponents for hydraulic radius and slope. Both the coefficient and the velocities © 
- are given as functions of bed forms and bed material. Velocities Bye by 


the formula check closely with those obtained in the laboratory. 


fet 
4 
Be nikal title of this research was “Analytical Study of Alluvial © 
Channel Roughness”, which was a research project granted to the first 
_by the National Science Foundation. The purpose of this research was to find ne 
a suitable formula to determine more accurately the mean velocity of flow, | 
and thereby the discharge of water, in alluvial channels. At the beginning of - Mg 
the research, the authors intended to study the variation of either Manning’ ae: 
Toughness coefficient or Chezy’s discharge coefficient as a function of the 
characteristics of the flow and properties of the sediment. It was found later 
that such an analytical approach is not likely to succeed. Anew velocity | i: 
formula was attempted, (1 the result of which is presented in this paper. 
Mie ing the mean velocity ofan 
pe stream more clearly, it is necessary to understand the mechanics of nf 
_ turbulent flow in pipes, in open channels and even in the turbulent boundary 
"layer, and also to understand the mechanics of sediment | transport. After 
"considerable review of literature the authors came to the conclusion that ree 
Note: Discussion open until April 1, 1960. To extend the closing date one month, a 2 
written request must be filed with the Executive Secretary, ASCE. Paper 2260 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
Society of Civil Engineers, ‘Vol. 85, No. HY 11, November, 
| resented at June 1958 Convention in Portland, 
1. Associate Civ. - Engr., , Civ. Eng. ‘Section, Co Colorado do State ‘Univ. » Fort rt Collins, — 


2. ‘Formerly Graduate Research Assistant of Civ. Eng. Dept., , Colorado Sta 2 


i Univ., . Fort Collins, Colo. . Graduate Research Assistant of the > Dept. of | 
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theoretical | approach to the problem cannot be obtained at ‘the present time a 
_ therefore an empirical approach was adopted. In order to facilitate the em-— 


correlation the method of dimensional was used so that 


| Literature Review _ 

Most literature concerning the equations of estimating mean velocity or 
_ discharge is for | clear flow in rigid conduits, , either pipes or open channels. = 

About 1768 Chezy(2) proposed a method of estimating the mean velocity of a 

by comparing: the flow conditions w with those of mitten i having 
Such a has been customarily written | in a form kn 


draulic radius, and S the slope of the channel. 


which a, b, constants and n n a facto 
a 1889, Manning(3) proposed several fo formulas for | estimating the mean 
velocity of turbulent flow in conduits. The following well known Manning’s — 


in which | M is an constant upon the 
_ However, Manning did not recommend its use because the equation is not di- m3 
tonaliy. | (3) is rently w written for the English 
which | n is 3 the Manning’ s roughness | factor. ot 
2 Re The authors introduce these commonly-used formulas here as reference | 
showing the variation of the exponent of the hydraulic radius and that of the id 
_ slope. Additional information regarding empirical velocity formulas can be 
: - found in the book “Hidraulik” written by Dr. s. Kolupaila (4) Kolupaila shows 


pends upon the velocity distribution, which is related to the mechanics ii 
turbulent flow along boundaries . The equation of motion for turbulent | flow is a 
known as the Reynolds equation which differs from the common form of the 
Navier- by additional called: the stresses. ‘The 


Bec 
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faa the three Reynolds equations together with the equation of continuity 
q for turbulent flow are not sufficient to determine the Reynolds stresses, ad- . 
md ditional equations must be obtained either through hypothesis or through ex- a 
perimental measurements to evaluate the unknowns. 
Among the various formulas of velocity distribution proposed for turbulent 
~ flow, the logarithmic law (see Eq. (5)) is frequently used by hydraulic engi- ee 
¢ neers. A brief review of this law may be helpful to understand its limitation _ 
of application. The logarithmic law (Eq. (5)) can be obtained either from _ add a 

_ Prandtl’s hypotheses(5) of mixing length by assuming that, near the wall, the Pies oe 
length is linearly proportional to the distance from the wall and 
= stress is constant, or from Karman’s similarity hypothesis(6) by as- - 

suming that the mixing length is only a function of the velocity distribution aw 7 

_ and the shear stress is constant. Therefore, this logarithmic law is for _ os ie 

‘turbulent flow n near It can be as\") 

in which u is the local mean velocity along the flow direction at a distance eS 


from the the boundary, V V, is the ‘shear ar velocity y — in which AT? is the local 


boundary ‘shear, is the so- Karman universal constant and the value 

of Z, is dependent upon a length parameter indicative of the hydraulic con- tt 

dition of the boundary. 

a From Nikuradse’s data for turbulent flow in pipes, it can be found that in 
i case that Vk/v is less | than about 3.5, in which k, is the size of the sand 
used in the experiments, the boundary can be classified as hydraulically 7) <a 
smooth and Eq. (5) can be written for flow outside the laminar sublayer — 


es ae found that the logarithmic law is not applicable to the flow near 
“the center of the pipe, which is self evident according to the assumptions used ; 


in the derivation of the law. If the logarithmic law were exact to describe = 


velocity distribution of turbulent flow in pipes, the total discharge, and hence 
_ the mean velocity of the flow, could be determined by integration through ie 
“use of the logarithmic law. It was found that the constants in the, resultant = 
"equations have to be modified in order to yield satisfactory results. ci at + 
In general the formula of mean velocity for turbulent flow in a smooth pipe Pee 
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_ Nikuradse formula for pipe flow can be applied to open channel ‘flow. How- ; a 
ever, Powell(10) found that because of the existence of a free surface in the “53 
open channel flow such an extension of Nikuradse’s work to open channels can- 

not be done successfully. Additional information for flow in open channels 

composed of artificial roughness. element on the boundary can be found 
works of Albertson and Robinson, (11) Sayre(12) and Johnson. 

.. the e foregoing r review of the logarithmic. law, there are two points which 
are important to the present study; (1) : although the. logarithmic law for turbu- 


lent flow near rigid boundaries has been verified by experimentation, Ps 


'y roughness is in accordance with apt 
Millikan(14) raised some doubts about the Karman-Prandtl hypotheses and 
mreed that without employing these hypotheses the velocity distribution of a 
— flow in pipes or channels follows the logarithmic law in the overlap — 
"zone where the “law of wall” ¢ and the “velocity - -defect law” are both applicable. 
the law of wall, , which is due | to Prandtl by use of a dimensional analysis, c 
is essentially first enunciated in 


general form by Karman can be written as 


= Umax at z and] h is the value o of at pipe or a. 
t 
discussion on will be later in con- 
- nection with the review of the turbulent boundary layer. The following re- | 
marks may be related to the classification of boundary roughness: nila 
“Shae the case of a rough boundary, the effect of viscosity on the velocity Bis 

can be discharge formula, such as Manning’ s which 


the case of. a the effect of the roughness elements 
e velocity distribution can be maglacted:: In addition to Eqs. (6) and (8), nd 


the -power velocity-distribution law: 
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a and D, i is the diameter of the pipe, and f is the Darcy- -Weisbach resistance ae 

oefficient. . According t to Schlichting, (15) the exponents in Eqs. (12) and (13) 

al are not constants as the Reynolds number increases, but dependent upon the. ie 
" Reynolds number of the mean flow. Eq. (13) can be written as an exponential _ 
_ type of discharge formula which will be discussed later. 

Since the problem of the mean velocity and the velocity distribution of 
a, turbulent flow in open channels is essentially one of a turbulent boundary © eo 
? ae a brief review of the literature on the turbulent boundary layer along a 


flat plate at constant pressure n may shed | some light « on the } problem of velocity 


a From extensive wind tunnel measurements it is found that the mixing length ei 
theory has many limitations and inconsistencies. At the present time, scien-— 
tists seem to be in favor of using statistical mechanics to study turbulence. 
A - completely satisfactory theory of turbulence is not available; scientists are 
eeking for laboratory data so that some satisfactory theory of turbulence can 


; ‘-, It has been found(17) in the wind tunnel that the logarithmic law is valid “pv 
only within about 15 per cent of the thickness of the turbulent boundary layer. q 


"According to Clauser, (17) the flow within the turbulent boundary layer can be 
_ divided into two regions. Ih the inner region, the law of wall is applicable; in 
“the outer region, the velocity defect law is applicable. ‘In the overlapping — 

rs _ zone where both the law of wall and the velocity-defect law are applicable, ais f 

the logarithmic velocity distribution prevails, which is similar to Millikan’s 

conclusion for turbulent flow in pipes and inchannels. __ 

_ Also according to Clauser the inner portion of the responds to the 
_ wall shear much faster than the outer portion. While the inner portion com-= tis 

pletes its response within few boundary layer thicknesses traveled, the 

outer portion takes a distance of tens or even hundreds times of the boundary , 

7 layer thickness for a corresponding response. A comparison of the response — ie 

- distance and mode of response to disturbances of various kinds and intensi- _ = 

_ ties confirms that a boundary layer is a truly non-linear phenomenon. Conse- -— 
quently, progress cannot be made by applying a linear concept of predetermi- — 

= response distances or times. Since the outer portion does not respond 


to the wall shear very quickly, the velocity distribution in the outer portion ‘a 
_ depends also on the history of the flow. _ Although the law of wall has been i 
: found to be independent of the pressure gradient along the boundary, it has not 7" 
been proven to be applicable to the case where the boundary is movable or 
flexible, such as the case of alluvial boundaries. —__ irre aa 
ve In brief it can be stated that, at present, there is no satisfactory theory of | a 
turbulent flow available so that the complete velocity distribution in turbulent ae 
flow can be calculated or predicted. _ Furthermore, since the flow is non- AE d 
linear in nature, it is very doubtful that a theoretical and exact solution of the : 
turbulent flow problem will soon be available, even though an seRrermnte < 
solution may be possible after extensive experimentation. me 
Another important factor in the study o of the mean velocity in alluvial Pe i 
channels is the presence of appreciable sediment transport which is absent E. 
“from flow in rigid channels. For flow transporting sediment, there are 


involved: a: (1) the amount of sediment transport, (2) the 
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problem of roughness and its effect on the discharge of the 
_ ‘There are several formulas ((18) , (19) and (20)) for estimating the camel: of 
_ sediment transported On the other hand, there is very little literature pro: x 
a velocity formula for alluvial streams. In case the bottom is plane, the 
alluvial boundary has been treated as a rigid one. _ For example, ‘Strickler(21) 


proposed that the ‘Manning’ s roughness factor can ue expressed as a we * 


One o of the major difficulties in determining of 
streams is that the bed configuration usually changes with the flow, condition 
- Consequently the bed roughness, which affects the velocity, changes with the’ | 
te _ flow condition. In 19% 1950, Einstein and Barbarossa(22) proposed that the bounda- 
- ry ‘shear of a dune bed be divided into two portions: (a) that pertaining to the J 
# — grain roughness and, (b) that pertaining to the dune roughness. Although such { 


“ an a roach seems lo ical, its a lication to ractical roblems is still ver 


» shown that cause a re 

auction in the resistance claim that the discharge and sedi- 
_ ment load cannot be expressed as unique functions of the depth, slope and sand ‘ 
size. _ This point of view is shared by some investigators of the problems of | 
4 flow transition due to the sudden change of the bed roughness. 
fe In conclusion it can be stated that the theory of turbulent flow is still not 
_ Complete even for flow near rigid boundaries. Its development for the case 

oe of a flow near a | movable ‘boundary seems € even more r remote. _ Moreover, the 
a effect of sediment transport on the resistance coefficient is still unknown, 

an - Hence no theoretical analysis can be made at the present time, for the _ ae 
problem o of mean velocity of alluvial streams. Therefore corre- 


ation seems to be desirable for | engineering purposes. _ 


= _ Analysis of the Problem of Mean Velocity in Straight Alluvial Channels 
: — the case of flow carrying sediment, the change of flow causes not only " 
3 the change of sediment transport but also the change of bed configuration. ahioe 
‘The phenomenon of sediment transport can be described by assuming that: _ 
 & the bed material is granular and cohesionless, (2) the amount of supply of 
_ the sediment is equal to the amount of sediment transport, and (3) the flow is" 
urbulent, steady and uniform. Let the bed be initially plane at a small — 
harge with no sediment moving, as the discharge oars’ following — 


< 


e 
is further increased, ‘the casvelesat of sediment b becomes more re intense. 
It can be stated that statistically there is | a critical condition under = 
which the movement of sediment begins. 
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Plane bed—Movement of sand grains occurs by rolling and sliding- 
for 


Ripples—As the is increased s still further ripples appear on 


the | bed at a certain stage. A ripple bed is characterized by a rather 
i ae wave pattern. The amplitudes of the ripples are usually small | 


compared to their wave lengths. ‘The characteristics of ripples are 
(c) Dunes— —At a later stage, dunes appear on ae bed. a“ dune bed is nso = 
; _ characterized by a long upstream slope with a steep downstream slope. — i 


an change of bed roughness. 


j -On further increase of the discharge the dune is con- 
siderably modified. The dune-spacing is elongated, and the dune- crest, 
6 is flattened. - Some of the dunes are washed out. This bed form may be =f 

called bars. The roughness of bars is not as high as that of 


) Flat smooth bed— -AS the discharge is further increased a stage is ha i 
oa hi: reached at which the bed becomes approximately plane. The bed — 
«ness at this state is less than that of dunes. This isa transition stage fied 
between dunes and antidunes. Vanoni and Brooks(23) call this type of me 
rot form as flat bed. A flat bed is always formed by the e hydraulic 
ee action of the flow, its deposition is generally very firm. . For given bed ~ 
‘material the hydraulic roughness of a flat bed may be smaller than - 
ofa bed. Hence it is come flat smooth bed in this study. A flat 
(e) . Antidunes—In the antidane is the up- 
stream slope of a sand wave and eroded from the downstream face of | 
_ the sand wave, consequently the sand wave moves upstream while ‘ae 
__ Sediment is transported downstream, 
_ The various bed configurations can be estimated from Fig . 1(25) for a given fie 
flow depth, slope, bed material size, and fluid temperature. Fig. lis ‘copied 
from an article entitled “Discussion on Mechanics of Sediment Ripple For- _ 
mation”(25) by the first author. A thorough understanding of the development — % 


locity ¢ of flow. _ The two parameters — and -— in Fig. 1 are substitutes for — 
— and which were proposed originally. (27) This is s permissible. 
upon which the criteria for various bed configuration are based can be under- i 
stood more clearly. However, the shear velocity Vv, is contained in both the 
‘ordinate - ag and the abscissa ~*°, therefore a method of trial and error is re-_ 
quired in order to solve for it. _ On the other hand, if the pair of parameters 
and is used as shown in Fig. the mechanics upon which the criteria 


for various are based cannot be visualized no 
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7 ‘Shield’ (26) cx beginning of motion and Liu’s(27) criterion 


A a the beginning of ripples are shown respectively. A modified criterion for  be~ 


ginning of ripples and criteria for the formation of dunes and for transition — 
 wegedes by Albertson, Simons and Richardson(28) are also shown. It should 
:  - pointed out that Fig. 1 is based ‘mainly \ upon laboratory data fo for which th the 
depth of flow is usually small. In order to maintain a laboratory flow with | a 
; ee _ bottom shear velocity comparable to that of a natural stream, the Froude 
Y ‘ - number has to reach values considerably higher than that of natural streams. 
‘Furthermore, the relative grain roughness ofa laboratory flow is generally 
4 than that of a large natural river. Such factors are not represented 
in ‘Fig. 1. ~ The authors also found that for 1. lack of exact and unified definiti no 
a various bed configurations, data from various sources do not always as agree 
_ with the classification shown in Fig. 1. Fortunately most data used by the © 
_ authors in the analysis contain information on bed forms. In case the infor-— 
a. mation on bed forms was missing, Fig. 1 was used as a guide in sonteatie. 
yee ‘change of bed configuration is illustrated by Fig. 2. ‘The . curves are for 
_ pipe flow after Nikuradse, and the points are data of No. 8 sand pertaining to 
the movable bed taken from Report No. 17 by the U. S. Waterways Experi- if 
3 mental Station. n.(29) As long as the bed is s plane, , the variation of the resistance 
ie coefficient as the flow changes is similar to that of pipe flow. The sudden in- 
crease of the resistance indicated by the points occurs as sand 
ay waves appear on the bed. From this it is easy to see that application of the 
Nikuradse approach to the case of alluvial channels i is not to o be 
‘The exponential type of discharge formula is already in existance a flow 


_ in rigid channels. For example, the Blasius equation which is for ms 
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(water) 10-5 ft2 /sec. 
- Note that cia Eq. (16), a variation of temperature of 20° F from 65° F 
changes the mean velocity about 4 per cent. Hence the effect of 
a eel on mean velocity is normally negligible if the correct formula ae 
used. The Manning formula (Eq. (4) which is! for turbulent flc flow near a ‘rough — tog 
haat is also an exponential formula. Note that these two formulas, Eqs. ahs 
(4) and (16), are for extreme cases, and the exponents are not the same. It is sae 
- possible t that the e exponents « of | the discharge > formula fo for turbulent flow in the a oa 
transition region have other values. In general the « type of of 


“tower 


in which cr an em nical coe! icient, x are pure 
y Pp 


2... _ For two- siesta steady, uniform flow, the depth of flow D can be con 
pre dependent mainly upon the following variables: q the unit discharge. a 
of the flow, S the slope of the channel, p the fluid Gensity, u the fluid viscosity, — aM 
g the gravitational constant, d the 50 per cent size, i.e., the mean size of the — : ‘os 

bed material, ¢ the standard deviation of the size of the bed material, AYs the 
difference in specific weight between the bed material and the fluid, and 
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. a the effect of sediment mixture on the flow depth is considered t to Tale > secon- 
; dary, the term a/D can be omitted from the equation. 3 It is difficult to use 


5 


> (21) because there are four terms containing the unknown mean velocity 
to this ‘difficulty the transformation can b be made: 


dailies a new set of dimensionless terms can be substituted in Eq. ( 


? 


), pV, can be written as Tp the boundary shear at the bed vale) 
Fora | wide channel Tp is equal to YDS. For a narrow channel the wall-effect ] 
may be appreciable. Th th then becomes YRpS v where Rp is is the hydraulic 


where Tp = and can also” 
be written as — 
In the grain of the bed material is spherical th value of fo 


- —_ nage Since q = DV, and the fall velocity of the sediment particle w is a function | bay; 
as 
— of th 
4 | ing t 
one | 
— 
0.008 
| mate 
were 
were 
w= 
= such The 1 
aral 
| 
plac 
oF th 
— 
18 pel 
and moving in Water depends upon Only the grain size and the water flow 
jj viscosity. The term >- cause 
the following analysis, Eq. is used as a guide in the correlation of 
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Empirical Correlation of Data 
No specific for 
Existing data on water. transporting natural sediment were collected as much 
as possible. Both laboratory data and canal data were used. Although most 
of the data were for flow having bed load only, considerable data for flow hav- ; 
“ing both suspended load and bed load were used. The depth of flow ranged 
from a few inches to several feet. The velocity of flow varied from less than 
“one foot per second to six feet per second. The slope of flow varied from 
0.0004 to 0.028. The size of sediment varied from 01 mm to 80 mm (approxi 
mately 3-1/2 inchea). Both uniform bed material and graded bed material 
were used. The variation of flow viscosity and sediment density of these data 
were not appreciable. The effect of side wall was (wan coer according to 
As explained earlier a theoretical treatment of the problem of mean ve- .. 
locity is not possible at the present time. ‘The result of this study, therefore, 
has been obtained from empirical correlation based upon physical reasoning, es 
dimensional analysis, and the mechanics of boundary layer. The drawback of ae 
using empirical correlation is that the cannot be 


for flow ina narrow Channel where the wall-effect is appreciable. 


parameters, with the exception of Fr, are included in Eq. (23b). . Note that the — t 


parameter Fr can be written in terms of and § according to Eq. (22a). 


Therefore, only those dimensionless ot Eq. (23b) are 
There remains a question whether the flow Din _—— canb re- 
plac d by Rp. The ratio (or r v=) is inns as either the Froude number 
or the Kineticity of flow. It is defined as either the flow velocity divided by 


| the celerity of the gravitational wave, or twice the kinetic energy divided by © 


the potential energy of the flow. Obviously both definitions cannot be apetied a 

| to the However, the ‘Substitu on of Rb for L D Eq. 


is permissible, plots that the is used in connection with t the 


flow resistance. Then == can be as index for the loss 


V+ and th the energy the parameter 
= 
£ 

equivalent y /S can be considered as an index indicating the ‘energy loss” 
caused by the surface waves. 
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: ‘The first parameter is the shear-velocity Reynolds | number and the second 
_ parameter is the fall-velocity Reynolds number of the grain. The third Het 7 
can be abbreviated as K,i.e., 
may need some explanation. These dimensionless parameters were evolved 
7 from a pot made od the first author(31) . in his previous study of the roughness 
alluvial 1 beds. this earlier study only Ved were u used. 
‘The t interpreted as the tractive force divided by the sub- 


ght of the particle multiplied by a resistance coefficient v,/V. a 
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and were added to the 


2. To the data =" 


San Luis Volley Conne! 
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= the plane bed changes into wavy bed can be estimated according - a 
re Fig. 1 (ripples are considered as incipient dunes). The intercept ofthese 
a lines with the line of K = 1 depends upon the third variable wd/v - For 


_ from those pertaining to the dune bed. These s straight lines for a plane bed 
are. Fig. Their slope is 1:0.555 horizontal to vertical, and t the 


size of | 


size of bed material, in mm 


Sa 


with d 


4 _ any given constant value of wd/v, the data, when plotted according . 
|  V,yd/v against K shown in Fig. 3, were found to fall on two straight linesde- 
— 
>, 


intercepts w with K = shown ii in ‘Fig. 5 according to the third variable 
). A general cmenti can be written for 


which 2 is 555, of the third va able wd as 


-_ Fig. 5, and A, m and N are pure numbers taken from Fig. 6 which were ob-_ 
tained empirically as functions of the mean size of the bed material. It should 
be mentioned that the curves for A, m, and N have been chosen empirically — 
_s that data can be plotted on parallel straight lines shown on Figs. 4 and 7 mt 
_ in other words, if the values of A, m, n, are to be plotted on Fig. 6, they — e 


From Fig. 5, for > 1000, the factor A be 


in which ¢ has an approximate value of 0.39, and is dependent upon 1 the shape 
_ factor and also the fall-velocity Reynolds number of | the | arti of bed ma- 


hey 


is assumed to be spher cal, the parameter os is not an indepe dent 
ariable and can be omitted. ‘The inclusion of — in the data analysis seems 
‘This contradiction does t exist 


othe r ameters i in Ea, (24), 
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— Own mh Figs. 5 and 8, the coefficient A in Eq. (24) depends upon eet |, ae 
the basis of spherical grains. In the 

— 

14 used in computing 

the parameter wa Can be written aS a com | 

yen Eq. (25) is substituted in Eq. (24), the result 

shape is considered. When Eq. (25) is substituted in 


At high Reynolds number Cp is a constant, therefore, Eq. - (26) does not con-— my, 
“toe the factor of viscosity, u or v. Eq. (26) is for turbulent flow near a. 3 4 
boundaries. If the values of A, m, N, and Q for plane bed of very large bed = 
; a material, that is, Q = 0.555, m = 1/6, N = 0.6 -—_ ra = 0.2 are substituted ed i a 


Eq. - (28), the result will be ais 


uf 


n 


V = 3: 35 
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The exponents of d, Rp S are the same as those ¢ given n by the Manning 


0.39, Eq. (27) will become the ordinary Manning- -Strickler formula. Theo- 

4 retically both the coefficient « and the drag coefficient Cp are dependent = tae 

_ the shape and the fall- -velocity Reynolds number of the grain of the bed ma- ar 

terial. The fact that both « and Cp are of the same order of magnitude ae 
the fact that the exponent of is as ~ as that of Cp may explain why 

the mean velocity may be proportional to Cp* *25. The factorse andCpin 

Eq. (26) can be used to explain why e Strickler’s coefficient is different = ’ 


ii in which. disa constant, approximately of 3.4, p is the slope of the curve, ap- a — 


1/2 and w= 1/1888 for ‘spherical grains. Su  Substitut- 
g Eq. (29) in (24) with = 3.4 yields 


For the limiting values of A, N, m and 2 for plane bed composed of very i _ 


bed material, that is, Q = 0.555, N = 0. .8, m = = 1/7 and A= - 0.287, as used, the 


- "exponent of Rp reduces to 0.72, ‘the exponent of S to 0. 57, ‘the exponent of a i 
> jf reduces to zero at p = 0.492, and consequently the exponent of v reduces to aa . 
Eq. (30) is then reduced to the Blasuis equation for turbulent 

similar « correlation can be found for dune bed as ‘shown in Figs. 7 and 8. 

oO The data shown on Figs. 7 and 8 can be represented also by Eq. (24), except — 


for dune bed 4s 0. 565, and the ve and N for bed 


of alluvial streams. Although Eq. (24) is dimensionally 
homogeneous, it is not convenient to use. a further simplification of ™ cy 
is discussed in the next chapter. 
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or ‘Ea. can be reduced 
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te 


in which P is P is not th the used in | Eq. (29), but is 


d, Ps, &; and and Q, , which depends upon the bed Ca to 
p 
more, the bed material is generally composed of a mixture of non- spherical — 


grains and the coefficient C, depends also upon the shape factor and the rine a 
— deviation of the bed material. In Eq. (31) x and yi are pure numbers 


vets 


VARIATION OF X 
BED 


func 


= 2b) the 
— 
2c) and, 
a : 4 in which A, m and N are shown in Fig. 6 a WE size 

— 


"When the values of Q, 


shown in Figs. 9 and as of the bed 

the size of the bed material. Because A, m, and N are empirical, the curves _ 
for x and y are also empirical. Since x and y are pure numbers, it is reason- — ce 

able to assume that they depend upon some dimensionless parameter, rather or 

- than upon the size of the bed material alone. The dimensionless parameter — 

_ which is still 1 unknown should be directly related to the boundary | conditions, 
and/or to the flow conditions. _ The unknown dimensionless parameter may be a 

some combination of those given in Eq. (23), with the possible exception of a 

oe i, _ Such a dimensionless parameter has not been attempted by the authors. 

_ On the other hand, since the choice of x and y depends partly upon bed con- 


figuration which is governed by two dimensionless parameters, V Vx/w, yy wd/y 


and the later can be obtained by dividing ~$" by 4, the effect of the hydraulic en 


_ boundary condition on the choice of x and y has been partly, if not entirely, 


considered. Further research is needed to me exponents x and y as 
functions of certain ‘dimensionless parameters. 1, qi 
a _ Note that in Fig. 9, the variation of the qxonsat x against the bed material 
“size d for dune bed is opposite to that for plane bed. _ Both the exponent it x for 
"plane bed and that for dune bed are 2/3 when the bed material size d is 5 great-— 
er than 4 mm (the exponent x of 2/3 is the same as appeared in Manning’s 
ote ) For plane bed, : x increases asd decreases for d ‘smaller than 4 mm. 


ae Note that in Fig. 10, the variation of the exponent y against the sediment 
7 ri d is exponent y for for 
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an bed | is 12 when the bed material size is greater than - mm ae y-value 
; a of 1/2 is the same as appeared i in 1 Manning’s s formula). — As the size of the bed 4 
a material decreases, the exponent y for plane bed increases for d smaller than | 
4mm. The y-value reaches an upper limit of 0.57 whend becomes 0.1mm _ 
or smaller (the y- value of 0. 57 is the same as appeared in the Blasius formu- 
ih) For dune bed the exponent | y becomes 1/2 when the bed material size is 
20 mm or greater. For d smaller than 20 mm, y decreases as d desceahen.” 
; is y-value is 0.30 when the d-value is 0. 1 mm or smaller. Meee g 


ee ‘That the exponent y y for dune bed a * below 1/21 may need some dis- | 


soi 


For turbulent flow near rough boundaries, 6 = therefore S is proportional 


In case of 6 is than | zero, S 4, such as in the case of 


's regime theory, the mean velocity ina regime channel is pro- 


‘portional dient to the one-third 


G 


"ge 
— By examination of Fig. 9 and Fig. 10, jee concluded: the Manning 
by _ formula is applicable for plane bed when the size of the bed material is 4 mm 
aie or greater, (2) the Blasius formula is applicable for plane bed when the size ae 
te of the bed material is 0.1 mm or smaller, (3) when the size of the bed canes 
ve al is 20 mm or greater the Manning formula is applicable regardless of the _ 
: oe bed configuration, because the effect of dune formation, if any, on the mean aly 
% <4 velocity is negligible and (4) the formation of dunes generates additional cfs 
ate en gy loss so that the energy loss is proportional to the velocity with an >= 
The coefficient C, can be from Eq. (32). However, sucha 
method of determining C, is very tedious. _ Instead the average Ccmiieet | ' 
were determined by substituting available data in Eq. (31) through the use of 
wwe 9 and 10 for choosing xX and y. Since Cc. is not dimensionless, and its 
dimension depends u upon the deenend of the hydraulic - radius, therefore the fe 
-Ca-value for the English system (Fig. 11) differs from that for the metric . 
q system (Fig. 12). The coefficient Ca should be a function « of the properties o 
a a _ the sediment and the fluid. Under ordinary conditions of sand transported ri 


ye water flowing in open channels, the density of the sediment is approximately 
See _ The shape of the sediment particle and the properties of water can 


also be considered approximately ‘constant. Therefore, Ce is essentially a 3 
function of the sediment size alone. It should be noted that the temperature 
~ variation of the data was between 15° Cc and 30° C (59° F and 86° F respectiv 


ly). It was found that the variation of Cc. due to ‘temperature change is les 
than that due to error of measurements. That the effect of the variation o 


viscosity of the water on the mean velocity | is small has been shown to be true j 


in the Blasius formula. ‘Therefore, the effect of the 2 temperature on the 
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“acharge can be neglected for practical that 


the appropriate exponents for R, andSareused. = 
Note that the effect of temperature is included in Eq. (24), which has been aa 


4 t _ reduced to Eq. (31). However, the effect of temperature is neglected in the 


_ determination of C a-curves. _ The explanation for this may be that, if the - 
| viscosity factor is used to form a | certain dimensionless parameter in nthe — 


complete. On the other hand, for practical purposes, the effect of neglecting 


viscosity rate the Pere — may be small compared to errors of other 


pa Because the exponents for the hydraulic radius or the slope Sdepends 
upon the bed configuration, separate Ca- curves for plane bed and dune bed are 
also ‘necessary as shown both in Fig. 11 and in Fig. 12. Net eae Fe 
_ In the case of a plane bed, the discharge coefficient Cg decreases as the — a. 
bed material | size d increases except when | d is smaller than 0.1 mm, , then the 
Ca-value is essentially a constant equal to 287 (see Fig. 11). This is the case 
Benne the Blasius equation is applicable, and where the discharge coefficient Bs 
is of the height of the boundary roughness. Between d = 0.2 mm 
and = 1 mm, the value of Ca decreases rapidly 2 as the value of d increases. 
‘Between d = 4 mm and d = 80° mm, the value of Cy can as, 
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Eq. (37) to the shown as as ‘Eg. a6). Notice that 
the drag coefficient Cp of a sphere decreases abruptly at about — ary 


which corresponds approximately to 80 mm at 20° C. Therefore Eq. (37) 

and Figs. 11 and 12 should be applied with caution for d larger than 80 mm. ll 

ey val In the case of a dune bed, at the bed material size d equal to about 0.2 mm, 
ww discharge coefficient is a minimum, which may b be interpreted as that the 

effect of dunes on the discharge coefficient is the greatest : at about d = 0. 2 i 
ee mm; and at the size of less than 0.04 mm, the discharge coefficient is es- ir 

sentially constant at 21. The discharge coefficient Ca increases" as the size 


of the bed material increases from .2 mm to 7 mm, which means that within © 


i * = range of the bed material size the effect of dunes on the discharge coef- | 
ficient decreases although the size of the bed material increases. - The Ca aa 


curve for dune bed coincides with ‘that for bed at d which 


mn «It was pointed out earlier that the bed configuration after the beginning of 
motion can be classified as plane, ripples, dunes, bars, flat, and antiqunes. 


ie insufficient data. The most important concept out of this research is that the $ 


_ flow over an alluvial bed is divided into two classes: that is, the flow with a | 
cq ; plane bed, and the flow with a dune bed. The discharge coefficient C, and the 
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the stage of the dunes, the the exponents x and y for the 


bed have been used in determining the Ca-value. The Ca-curve for ripples" 
and sand bars is shown immediately above the Ca- curve for dunes. The two 
curves coincide with each other for the bed material sizes less than about 
0. 06 mm, and greater than 1.6 mm, which means the distinction between | 

_ ripples, bars, and dunes as far as their effect on the discharge is concerned — 
is nil for d< 0.06 mm and for d >1.6 mm (It can be estimated from Fig. < San 
that ripples will not form when d is greater than about 2 mm). ‘ That | the a 


- curve for ripples and bars lies above that for dunes means that the resistance 


re of ripples or sand bars is generally less than that of dunes. The minimum © e 7 = 


” In the case of a flat bed, which follows the stage of sand bars, the values a 
| of the exponents x and y for plane bed have been used in determining the Ca- 


Ca-value . for ripples and bars is about 18 at d equal to about 0.16 mm. ee oe 
_ value which is shown immediately below the Ca- -curve for plane bed. The Se F 
| value reaches essentially a constant of 256 (see Fig. 11) when the bed materi- cf 
al size is less than 0.1 mm. The Cq-curve of flat bed coincides with that of | 
plane bed when the bed material size is equal to or greater than 20 mm. That 
the C,-curve of flat bed in general lies below that of plane bed indicates a 
the discharge coefficient of a flat bed is normally less than that of a plane 
bed. ‘This means, for given bed material, depth of flow and slope of the ¥, 
channel, the discharge or peteay of a flow with flat bed is less than that of a 
a7. 7 It should be ashe that in order to select the exponents x and y and the iti 7 
propriate discharge coefficient C,, it is necessary to estimate the bed configu- 
-- ration, which requires the use of Fig. 1. _ The following is an illustration | to 
estimate the mean velocity of from Run No. No. 5 - 


0010: Report No. a, (29) mr 


‘Required: 


the sediment is w=0. 234 ina and 30.9 therefore 
cy = =0. 483. From Fig. 1 it is estimated that ‘the bed would be a dune-bed 


Therefore, the sa and coefficient for a dune-bed at d = = 0. 483 are ar 


| _In the case of ripple bed, which is considered as the transitional stage be- = (7 
Tb) | 
at a 
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_ The data computed in this manner are shown in Fig. 13 and Fig. 14. ‘In a 7 


_ el 13 which is for a plane bed, 83 per cent of the data are within 10 per cent 


oe of scatter. None of the computed velocity exceeds | 20 per cent of cnt om 


18 


| = eee the measured value. In Fig. 14, which is for dune bed, 70 per cent of 


a error is about 10 per cent or less. In view of the fact that it is very diffi- 
= to obtain accurate data of flow in alluvial channels, such as measuring vical 
depth of flow and the energy slope, an average discrepancy of 10 sont, 

locity 


1 % error line 
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4 ‘This 
facto 
case 
are. | A comparison of computed meon velocity with 
measured mean velocity for alluvial plone bed 


‘The results presented are for flow in straight channels. 
‘The effect of a side wall should be eliminated by using standard pro- | ve tt is 
cedures. (30,35) This involves a method of trial and error since Rp cannot be 
computed without knowing the mean v velocity. -Inor order to use this method, | ey 
bes the total discharge and the depth of flow must first be assumed. The 
mean velocity can be found according to the equation of continuity, then Ry ie 
can be ‘computed. By using Rp, the slope of the channel, and the mean size of 
_ the bed-material, the mean velocity can be checked according to Eq. (31). = a 
‘This method should be be used through re repeated | trials — the result is is satis- 


Eq. (31) is suitable for a steady, uniform flow. However, this is not a a 
case for most of the natural streams. . In the case of natural streams the dis- 


charge | coefficient has to be modified to suit the The 
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certain methods of empirical so doing it was 
wo sary to make some assumptions for simplification. However, in order to bee 
a the problem thoroughly so that the final solution of mean a 
> of alluvial streams | can be obtained, some additional research work definitely 
is needed. The needed research is almost unlimited. _ The following aad 


| The information on bed configuration is very important in order to ap- 


"ply the discharge formula properly, the classification of bed configu- 
ration needs to be defined more accurately, 4 
2. In this study the effect of sediment shape was not considered. All sedi- “| 
: ‘ment particles were assumed to be spherical. In order to improve the 


accuracy of the method, it is necessary to determine the effect 2 en 


3. The effect of the mixture (size gradation) has not been eiveatliiees 
thoroughly, further research is needed to determine its effect on the 


velocity. ~The accuracy of the authors’ method depends consider: 


upon the size of the bed- 


in Fig. 11 and Fig. 12, should be classified more "accurately ; according 


aac 
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1 ae Dynamic viscosity of the fluid 


fluid 


Kinematic viscosity of the 


Density of the fluid tp 


Standard deviation of the size of the bed material 


Local shear at the bed level, Tp = YDS for uniform flow in 
a very wide channel, and Tp = ’RpS for uniform flow in “hi 
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Dimensional exponent pertaining to Eq 


«Os 0. 555 for plane bed 


= 0.565 for dune bed 
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EFFECT OF AQUIFER ‘TURBULENCE ¢ WELL 


ieee change tines laminar flow to turbulent flow, in the he case of water vi 
ing through sand, occurs over a wide transition zone beginning with a Reynolds _ 
number of about 10. A method for ' estimating the exponent « of the velocity 7 
term for head loss calculations involving flow into wells is presented. _Typi- e. 
cal calculations show that head losses due to aquifer turbulence may have 


oa turbulent flow is present or may vt expected to occur in a given well. 
>’ Turbulent flow, as opposed to laminar flow, may be defined as flow under — = 
conditions that relate the energy losses to a velocity exponent that is oe il 
than one. In laminar flow, the energy losses are directly proportional to the 
velocity. Most well hydraulics formulas are based on laminar flow con- a 
"ditions, 4 i.e., doubling the drawdown in an artesian well doubles the yield. | par 
oe Many investigators have shown that the Reynolds number can be used in | 
"porous media to indicate turbulent flow. The usual definition of ans Reynolds 


4 anion request must be filed with the Executive Secretary, ASCE. Paper 2265 1 is ii si oa 
9 part of the copyrighted Journal of the Hydraulics Division, Proceedings of of the Ameri- a 
ean Society of | Civil Engineers, Vol. 85, No. HY 11, November, 1959. Bees oo a 

l. Field Engr, » Reward E. Johnson, Inc., St. Paul, Minn. 
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‘Thus the velocity is calculated by dividing the flow rate by the e gross area of 
formation through which this amount of water is passing. 


- sieves and N is the number of grains of diameter dg. N may be calculated - A i 
Been the grains are all spherical and have the same specific gravity; 2 bee 
.. thus, the number of grains per gram of sand can be calculated, and from this, = 6 

8 is obtained by multiplying the number of grains por gram by the number of 
zy ‘grams retained between the two consecutive sieves. As can be easily seen, — a fort 
ss The kinematic viscosity can easily be obtained from published tables if the 

Other investigators have indicated that if the Reynolds number is greater = 
than 10, turbulent flow rather than laminar flow exists. In some cases, two © D 
i (the equare) has been arbitrarily chosen as the exponent of the velocity after a 
_ the existence of turbulent flow has been established by calculating the ot i meal ini 
eee Reynolds number to be larger than 10. _ This means that doubling the flow rate Q 
causes four times as much drawdown head loss through zones o of 
a Tne calculation of the Reynolds number would be much simpler if so 
method of estimating the average grain diameter, other than by use of th 


Bie ss above formula, were available. _ A method, used by the organization with which 
ne iat the writer is connected, is to assume that the average grain size ofa water- | 


cent coarser. 


4 


scribe the sample as to coarseness oe was selected because it usually repre- 
fe sents the portion of the sieve analysis curve that has the steepest slope. Al- 
ee though this assumption is solely empirical in nature, the magnitude of the 
; ee Reynolds number is, in part, defined by its method of calculation and, there- — 
fore, (d) could be chosen by any reasonable method. The simpler empirical | 7 
% definition | of average grain size has been used in the investigation reported in 


To determine the Reynolds number at | which the head loss ceases fh ne 


Cumulative pe 


_ tain the exponent to be used for turbulent flow, three typical graded sands on - 
were subjected to laboratory tests. These sands are designated as A, B, and 7” 
Cc in the accompanying text and illustrations. Sieve analysis curves of the oe 


‘Sands are shown in Fig. 1. 1 In addition, data reported by Lockman, (1) pertain- J 
a to well-sorted sand and gravels that are normally t used in gravel packing 
of wells, and curves developed by Hudson(2) and Machis(3) were studied. 

Each of the tested sands was placed in a constant head permeameter wnere 
= he flow rate could be varied to obtain Reynolds numbers ranging from within 4 “ee 
tae) laminar flow range to the maximum number obtainable with the equipment.> _ 
Flow rate, head loss, and temperature were carefully measured. The © as 


Reynolds number was then calculated by using formula the average 


4 
© 
Dearing sand lormation is represented by the diameter Of sieve Upon Wich 
= 
; 
— 
qi 4 
— 
— 
— 
— 


from the sieve curve. The velocity was measured in and 


kinematic viscosity in ft2/sec. The addition of a constant for dimensional 
~homogeneity resulted in the equation: 


After the were vers tebeleted, a graph of the hydraulic vs. 


sands used in the experiments. 
‘Test of Sample 


e nolds number was constructed 0 logarithmic paper for each f tne three 
y' n log r ea i e 


The average grain diameter of this. poner as represented by the 70% re- 

ained size, was 0.055 inches. It is a typical coarse sand and fine gravel 


_ formation. The sample was tested under hydraulic gradients 1 ranging from a a 


100 


coar 


Grain size in thousandths inch 


pes denoted 


£ 
(2) — 
— 
— 
of 
= 
lis, 
rtain- 4 
king | 
by letters A,B,andC). — 


i "ranged from about 1 to 39. Fig. 2 shows the logarithmic graph of the aperenit 
— adient vs. the Reynolds number. The slope | of the first part of the curve is i 
unity. This indicates laminar flow, since the hydraulic gradient is proportion- 

et al to the first power of the e velocity. (Note that in flow through a a given sample 


of sand and gravel the Reynolds only wita changes in velocity.) 


033 inches. Water was through this material under gra- 


1 Bie - dients ranging from 0.4 to 7.0. In this experiment, the calculated Reynolds 4 
number ranged from 0.4 to 24.5. Fig. 3 shows a graph on which» 


This sample has an average grain di diameter (10% size) of 0.011 inches. . 2 

a | classified as a medium sand. The: sample was tested under hydraulic era- 
e about 6. A logarithmic. graph of the hydraulic gradient vs. the Reynolds mS 
a. number for two separate runs is shown in Fig. 4. This diagram shows two i 

lines, roughly parallel, one for the first run and one for the second run. a The t 

later data were judged to be more accurate as some air bubbles were observed 

in the test apparatus during the first series of readings. s. This w would tend to. > a 

increase the slope as indicated by the plotted data. The slope indicated by the §- 

curve of of test | run No. 2 is close to unity, indicating that the gradient was di- 5 
rectly proportional to the velocity throughout the test. ‘This result is to be 


expected since the maximum of the attained in this 


a 
The experiments conducted a Lockman(1) were made for the purpose of 3 
testing sand movement into gravel packs surrounding well casings. The data we 


and the results used to examine the effect of turbulent flow on head losses. = 
Sie (a) A very uniform sand, all passing the 40- mesh and all retained on the 


60-mesh sieve, was tested under four « different hydraulic gradients. 
These when plotted versus the Reynolds numbers on 


— 
— 
— 
f the curve is about 1.2 
4 ima initial straight line. The slope of the ws aulic gradient is proportional to the 4 
__which indicates that in this slotted perio fall slightly above the 
ani .2 power of the velocity. The last four Pp ld be greater at higher Reynolds — a ug 
—— | hich may indicate that the siope would be greater it ey 
— 
> 
itv. curves are es- 
t ht until a Reynolds i 
which the slopes become steeper (about 1.2) prem 
4 a = head loss) varies with the 1.2 power ofthe velocity. Sage 
— 
— 
— » 
— 
— 
— 
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showed a variation of the 2 gradient with the velocity to the 1.3 tonal 


The Reynolds numbers ranged from 6 to 46 in this series of tests. — 


A sample of similar material was tested in the laboratory maint comp! 
the s company. This material had an average grain diameter 
0% size) of 0.208 inches~comparable to the 0.203 inch material used = “indies 
by Lockman n. The hydraulic gradient appeared to vary with exponents” 7 | once 
ae = of the velocity ranging from 1.0 to about 1.3, through a range os ‘pened 

c) A 1/2- inch gravel tested under four hydraulic gradients y Lockman 

showed the possibility of two curves. The hydraulic gradient varied 

ve 
know 


trans 


with the 44 power of the ne arp ina range of Reynolds numbers from 


the friction versus the Reynolds number was plotted ye logarithmic in 
paper. The Reynolds number ranged from 0.1 to 100. Average grain sizes 
ti he calculated by formula | (2) thus the magnitude of the Reynolds numbers are 

slightly different from the values computed for the tests described above. 

i The curve obtained by Machis has a slope of unity until a Reynolds numberof — > 
is reac ed. Between Reynolds numbers of 10 and 100, the slope of 4 

curve changes. constantly indicating a transition to turbulent flow. Ata 
ie Reynolds: number of 100, the ae of the curve shows that the apiraniic 


Hudson(2) a graph showing the relationship between the fiche 
factor and the Reynolds number. This graph contains data from five separate a Tt 
_ investigations and, therefore, offers broader coverage than any ofthe curves f tuorby 
reviously « discussed. The range of Reynolds numbers is from about 0.1 to a ‘The | 
over 1,000. The initial slope of the curve is unity. Beginning with a value of ani 
the Reynolds number of about 4, the slope gradually decreases until a slope of 
about 0.4 is obtained at a Reynolds number of 100. a (This indicates that the eG of the 
gradient is proportional to the 1.6 power of the velocity.) me "grain 
_ Reynolds number of 1,000 the slope is about 0.2. This indicates that the re- na | analy 
lationship of head loss to velocity varies through a wide transition in zone with © : 
the exponent of oad 2 as a limit for of 


ns from Anz lysis of Test Data Reyn 
"When all all of the test data compiled in the of the company 10 
thea 
ae are compared to the results obtained by other investigators, ,it appears thata ff 
_ Reynolds number of 10 (as calculated in this paper) may correctly represent — _ veloc 
_ the limiting point between strictly laminar flow and the beginning of the tran- - great 
_ sition to turbulent flow. It furthermore appears that the hydraulic gradient be re 


varies with some power of the velocity that gradually increases from 1 to 2, 


as from 10 to infinity. 2 Hudson’s curve show 
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"that of the is still ata Reynolds num 


1 ,000—although the increase is at a lower rate. Therefore, it appears that the Te a 


- transition zone, from ‘Strictly laminar flow to the development of complete = F 
turbulence, is a wide one and that - rarely ev ever (at least in well problems) will i 

ee _ For practical engineering analysis | of well problems, the foregoing data oo a 
indicate that laminar flow in porous media occurs when physical conditions = 
are such that the Reynolds number is 10 or less. Under these conditions = 
head | loss varies directly with the velocity. If the Reynolds number is above cota, 

2 10, conditions are transitional between laminar flow and complete — 

| Here the head loss varies with the velocity raised to a power between 1 and 2. 
ep To simplify calculations and at the same time introduce a reasonable factor 

- of safety, a table was constructed for the purpose e of selecting the power of 

the velocity that the head loss may vary with when the Reynolds number is — 

x} 


between 80 and 90 


Conditions low ‘around and into a well the Reynolds number may 


The Calculation of the Reynolds Number and Its Use in Actual Field Practice _ 
‘a area immediately surrounding a pumping well is the critical zone for = . 
- turbulence because of the higher velocities of flow attained in that region. et 
The highest velocity is attained at the point where the water enters a well en 


screen, or where it enters the bore hole in the case of a well drilled in rock. - a 
- ‘This velocity m may be calculated by dividing the pumping rate by the total area a 
of the screen or the area of open hole exposed to the aquifer. The ome ar 
= diameter is taken from the cumulative percentage curve of of the sieve 


analysis at the point v where 10% of the sand, by weight, is coarser and 30% - 


of coarse sand similar to sample A.’ The yield of the ‘would 

have to be 370 gpm, in order to develop flow conditions corresponding to ee 
Reynolds number as great as 10 at the well face. . To attain a Reynolds number ~ 
of 100, the yield would have to be 3 ,700 gpm—a highly improbable figure for - =) 
the assumed conditions. The effect of turbulence would be represented by a ae 


Velocity term with an exponent slightly larger than unity for pumping rates 


greater than an 370 gpm, but almost certainly never as large as 2. Also, it alt 
be remembered that the increase in head loss due to turbulence occurs only =a 
a short radial distance from the well. y if the well in the above ee 
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at: distance of face it would be less 10. Fig. 
ince shows a graphical picture of a section of the cone of depression for this well. 
The graph is drawn on semi-logarithmic paper since the head loss (or draw- 


down) ) is proportional to the of the of the horizontal 


a eon well— 11 inches from the face of the 12-inch well. The cone of omesere 
‘from that point inward is distorted to the extent shown by the difference be- moi 
tween the laminar-flow dashed line and the turbulent-flow solid line. The in- D 
creased head loss due to turbulence is only 1. 4 feet, based on the use of an thy 
_ exponent of 1.2 for the velocity. _ The extra head loss due to turbulence repre- 
sents an increase of about 2% | over the drawdown calculated on the basis of | 


Be _ This example assumes no increase in permeability resulting from develop- ' 
i E ment. Inasmuch as there is some development in almost every well, its effect 
__ should be considered. As an ‘example, assume that the permeability within 11 

_ inches of the well face is increased three times by development. The average 
“s grain size in the developed zone will also be increased. It will almost double | 
at the well face and the resulting Reynolds ‘number will be about 54 instead of 

27. However, because this is less than 100, the head loss is still proportional 

. to an exponent of the velocity which is less than 2. 0. Assuming that it is 1.5, 

; bs the dotted line in Fig. 5 represents the profile of the actual cone of depression 

after both the effects of higher permeability and higher Reynolds numbers — 


F a a been accounted for in the developed zone. _ The total incremental draw- 


ing 


down in the zone of turbulence is 3.4 feet. Of this, 1.1 feet is attributable to 
the turbulent flow effects as 2.3 feet would have resulted from laminar co 
oe throughout the zone. . When this 1.1 feet which would be caused by the total. 
‘ mg drawdown outside the well is compared to the total drawdown outside the well — 
ee (66.9 feet), it is seen to be an increase of only 1.6 per cent. It is obvious that 
in this well, turbulence, even | though it exists, effect on x. 


haa Because Sample B is finer grained than ila and ‘thus has a lower 
ae Reynolds number at any corresponding pumping rate due to its lower average 

: Fa grain diameter, a yield of 620. gpm would have to be obtained before Ges. ett 
_ Reynolds number exceeds 10 at the well face. Assuming a well with the same 

_ physical conditions a as Sample A, the effects of turbulent flow would be even _ 


a well with the same physical dimensions were drilled in a formation 


of the well. if it were possible to obtain 1,700 gpm from a 12-inch well in ye 
this material, the drawdown would be in the neighborhood of 700 feet, an 
ra In coarse materials, such as the 2 1/4- -inch gravel ‘discussed in the report by | 4 
fuxbalence and the Reynolds number are more significant. At 


eynolds number would be 100 at the well face and less than 10 ata atmos 
from th face. 2. Fig. 6 as 


i a 
{ 
chart, the curve representing the cone of depression is a straight line. 
— 
— 
— 
4 q 
4 
— 
4 
a 
ite 
— 
_ reach a Reynolds number of 10 at the well face. Evenifa4-inch wellwere > | | ; 
constructed, the yield would be 400 gpm. Thus, it is unlikely that turtinlence 
— . 
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depress The average grain diameter is 0.208 When the effects of 
turbulence are considered without development, the drawdown is increased Bo 
a 62 feet above what conditions would be with | laminar flow. id ‘This | represents - 


this case, it is estimated that the average head loss will vary 
the second power of the velocity. _ This power was chosen because it repre- is 
Das the maximum loss conditions at the well face where the Reynolds number i 
is equal to 100. Actually the average head loss would vary with some ~~ i 
_ power of the velocity because the exponent of the velocity is increasing from | 
— 1.0 to 2.0 in the transition zone. . By choosing the exponent from the Reynolds © 
number at the well face a second safety factor is automatically incorporated. _ ; 
ae Suppose, on the other hand, that development reduces the head loss in the . 
gone of turbulence by 25% and that the Reynolds number is not increased be-— 
cause the average grain diameter is increased only slightly. In our case, the 
‘drawdown is only increased 0.47 feet over that which would occur under lami- 


- nar flow conditions and represents about 13% « of the total drawdown in the 


= 
— 


» or 

att! can easily be determined by the » Reynolds number criteria described » 

this paper. If the Reynolds number is 10 or less, the head loss ey with — 
first power of the the velocity and laminar - -flow well formulas apply. If the 

a calculated Reynolds number is between 10 and 100, the exponent of the velocity — 
can be selected between 1.0 and 2.0 respectively. It is suggested that an ex-_ 


ponent be selected | by determining t the Reynolds number at the well face (high- 
rofile 
i 


in a formation consisting of -inch clean gravel 85, =) 
Reret the contributions of turbulence to the t total drawdown were only 13%. By 
Therefore, it seems unlikely that turbulence in the formation adjacent to ~ = ) 
_ bore hole will alter the drawdown-yield relationships in most water wells. | < 
An important factor that t appears: to reduce the effect of turbulence is ‘maxi 4 

mum development of the formation materials adjacent to the well. It has been 

_ ttown that suitable development of a well in a 1/4-inch gravel formation could 


| 


5. Profile of the | 


the upward movement of water through the well casing. % All of these losses oi, 

_ probably vary with the square of the water velocity. _ Because the total losses: , 

are actually present, as shown from many field pumping tests and this paper — | 
indicates that only minor turbulent losses are present in the aquifer immedi-— 

_ ately surrounding the well, it appears that some of the above losses are more 


ers as it turns upward inside the screen; and the friction losses caused by " @ : 
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Other ground water investigators have indicated great losses attributable 
a we 

of f 
4 
— 
ii 


ae be readily calculated. a They are, for the most part, minor losses because > 


The screen losses are more difficult to calculate. However, it appears that 
‘the per cent of open area in the well screen is more important to high yield- — 


a 2. wells than previously recognized. If the well screen has insufficient open — 
area, not only are high friction losses caused as the water passes through the — 


screen, but even more pronounced may be the losses within a short horizontal 


distance from the well face. This is because the water must converge to the % 
limited number of openings in the screen and thus the gross area of flow with- 
in the formation is reduced which would make the calculated velocity and 
‘esulting number than what actually exists. 


- Lockman, John R. and Rohwer, Carl. er 
’ ‘Wells in Fine, Uniform, Unconsolidated Aquifers” . Unpublished Progress 
_ Report, Soil and Water Conservation Research Branch, Agricultural Re- — 


"Search Service, v. . Department of ‘Agriculture, Fort Collins, Colorado, 4 
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FORCES ON SUBMERGED STRUCTURES? 


binge 


F. BRATER,) ASCE, 


LESLIE D. STAIR.3—The authors wish to express the on 
“wae press their appreciation for the 
Belch of the discussors. The value of the paper has been ere 
a Mr. Alterman presented the results of an unusual field investigation which 
indicated that nara on piling can be accurately determined by means of the 
basic equations. He also presented formulas for predicting tt the force « on a 
horizontal cylinder. However, these are for a different type of motion than a 
_ that discussed in the paper. In accordance with Lamb’s presentation, (32) the 
~ force is caused by a uniform flow past a cylinder whose radius is ‘Small com- 

ared to its submergence. ‘The resistance is calculable from the rate at cone. 6a 
thich energy is supplied in the form of a surface wave downstream from the = a 

cylinder. — The derivation | holds for very deep water, and does not include t the — es 
c special case of the cylinder resting on the bed which Mr. Alterman presented. Bee 
| Furthermore, as stated in the reference, the force becomes zero in the limit- 
4 ing case discussed by Mr. Alterman for which C= Wed ‘gd, (d is the depth a 
_ water, assumed in the derivation to be very large.) The result is related m:.. 
the authors’ study only insofar as the motion within a wave approximates a 
if ‘Steady flow, and only f for that portion of the resistance which is caused by th the 
creation of a secondary : system of surface waves. 
- Mr. Hamlin called attention to the favorable agreement of the experimental 
mass coefficients with the theoretical ones” and with those determined 
Yu 19) for the barge-like models. ‘The question | of correspondence of c ef- os 
g m q on of co po Cc co 
ficients and the comment with regard to the frequency of the oscillations in 
_ the tests by yull9) ; are related to information which has appeared in several — — 
_ of the references in the paper,(5,6) in the discussion by Sarpkaya, and in _ : 
i more recent summary of the effects of unsteadiness by McNown and —it«*s 
woe Keulegan. (30) Mr. Hamlin is indeed correct in pointing out that the periods in 
; | Yu’ Ss experiments v were longer | rather than shorter (as stated in the paper) i 
Shmwnes those in the Michigan tests. The statements concerning frequency should 
have made use of relative frequency with the reference being the frequency of 


| vortex formation. the amplitudes in the by Yu were 


a. Proc. ‘Paper 1833, November, 1958, by Ernest F. _ Brater, John 8. McNown — 
and Leslie D. Stair. 
Prof. of Hyd. Eng., Univ. of Mich., Ann Arbor, Mic eet 
. Dean, School of Engineering and Architecture, Univ. of Kanens, | Lawrence, 
Kans. (formerly Prof. of Eng. Mech., Univ. of Mich., Ann Arbor, Mich. es 
Res. Engr., Boeing Airplane Corp., Seattle, ‘Wan. Assoc. Res. 
Engr., Eng. Res. Inst., Ann n Arbor, Mich.) 


Ba. 
ieD. 
3 
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riods. The 


in the tests of barges anc 
formation definitely played a role. 
-. Mr . Sarpkaya contributed interesting comments as well as new information. 
_ ‘The transverse waves could probably have been avoided by the method he ‘sug- ; 
gested. it was simpler to avoid the occurrence than a 
“ta Definitions of virtual and added mass can be presented or defined igs” 
several points of view all leading to the same result. The kinetic energy of 
the flow— ~and its time derivative—are customarily used, and the resulting 
‘derivation is entirely logical. . Darwin’ s concept of displacement is an intrigu- 
ing addition. Another form has been presented by Lamb, (37) in which the 
| integral of the Euler equation, including | the term for unsteady flow, is used. 


_ mental oe = is the product of an elementary mass and the local acceler- 


tis the any and U the velocity of the ol the 
Mr. Sarpkaya’s extensive results for separated and composite | 
te additional information which i is useful. His deductions concerning tests j 


are ‘seca valid although subject to some uncertainty, five he suggests. He q 


4 = has correctly pointed out the error in si 9 in going from Eqs. (4) and (5) to 

‘Mr. Herbich has as reported on 1 difficulties ‘encountered in wave 

tests in tanks because of the occurrence of transverse waves similar to those 


, - observed by the authors. Perhaps the suggestions made by Mr. Sarpkaya in 
- this regard and his references to the work at Grenoble will be helpful. __ 


0. “Vortex Formation and Resistance in Periodic Motion,” McNown, J. S., ye 


we ou 


age Keulegan, G. H., Jr. Eng. Mech., Proc. ASCE, Jan. 1959). bm, afne! 


37. "Lamb, H., loc ‘cit, 77. Yam 
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“HAROLD ZOLLER,! F. ASCE and ARNO 1 T. LENZ, 2,2 F. ASCE. —The 
ae wish to express their thanks to Messrs. Wright, Riggs, _ Schaefer _ 


; ie The work reported by the subject } paper was related to the Big Eau same 
* Basin because private funds were allocated for the study in this specific basin. 
‘It is true that many basins could have been found for which more precise __ 
‘meteorological data were available. interesting circumstance here is 
that plausible or at least explainable results were forthcoming even though 
the available data were transferred beyond normal distances. - Prior research 
on this basin(1) has predicted the total runoff volume to be anticipated during Ff 
any year as a function of snow cover , water table elevation and other ante- . : 


cedent conditions. _ Therefore, the question of quantitative runoff raised by 


_ Mr. Riggs was ‘considered by the authors although not reported in the subject = 


Ae The “additional variable” required by Mr. Riggs to explain and account for . 


the wide divergence of curves in Fig. 8, is the total anticipated runoff. ‘Figs. 
9 and 10 of the subject paper are intended to translate the total runoff ordi- 
nate for each year into a common base at a melt potential of 7 inches, becau: 
the runoff for each year wa was substantially complete at this melt potential for 7: 
Mr. Wright reported in his discussion, the limitation of the techniques 
+ | | ot as a means of forecasting daily discharge rate is the accuracy and a — 
| availability of such weather forecast data as air temperature, relative humidi- eo ‘ 
‘s ty, wind velocity and cloud cover. Only insofar as these data can be forecast 
_ reliably is it possible to make any reasonable forecast of discharge. a 
_ Mr. Wright indicates that the runoff forecast can only be made for one day in 
—* based on his study of the records for 3 years. The authors’ paper : 
made no attempt to forecast runoff, but it would be helpful if a forecasting _ 
technique might evolve from analyses suchasthis.§ 
oe i authors agree with Messrs. , Riggs and Wright that the method is only — 
_ applicable about 50% of the time me as a result of negligible s snow pack, frozen 
gages, or excessive precipitation during the snowmelt period for some years. = 
= discussions of Messrs. Riggs and Schaefer were concerned primarily % 


with the adjustments of runoff data for ice effect. Since the bulk of the ‘spring 


gees runoff usually occurs ina partes of from one to 5 days in the Big x 


q 

“ ‘ 

— 

— 

Prof. and chmn. of Civ. Eng., Univ. of Wis., Madison, Wis. 


-Pleine Basin, the authors were particularly concerned with the terminal ¢ 


= 6 revision is almost entirely a matter of pone Pecan since there are 
a no observations available regarding ice floes, ice dams, and the like. In his” 
discussion, Mr. Schaefer points out back water effects of ice on some 
* Wisconsin rivers as high as 5.69 feet. . The subject paper has not contended 
ay that back water effects do not occur, but that they are not likely to occur on a 
- the Big Eau Pleine River near Stratford because of the relative absence of 
channel obstructions. is interesting to note that the reference papers indi-_ 


cated by Mr. Riggs do. not include any ice index x values in excess of 3 feet and a 


subject paper. 4 


 Fuhriman ,D _K. , “Forecasting ‘Snowmelt Runoff” Thesis, Univ. of Wis., 
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ON SELF- -AERATED FLOW IN | 


ad torent 6. st G. Straub Alin 


LORENZ as and ALVIN 2M M. ASCE. 
In view of the well- known research on the subject of air entrainment carried : 
on at the University of Naples, the discussion of Mr. Viparelli in which he = 
compared some of their work with that of the present paper will aid ina ool 
under standing of the problem. There are many facets of the 2 air entrainment 2: 
_ phenomena in open channels that are still in need of investigation, —— 
the basic question of the exact mechanism by which air is insufflated into the 
flow. ms It was hoped to present ir in the present pager. data that would b be a — 


air were based | upon meamirements ¢ of for various 
flows and slopes. ‘The mean depth d and’ transition depth dy were defined in “Y 

| terms of the air concentration distribution. Knowing the water discharge and 
the mean depth which implies the area per unit width of water flow, the mean 

velocity may be calculated as was shown in Table II. ‘Sty 

_ The concept of the transition depth dy defined from the shape of the 

(concentration distribution curve provided | a means of establishing the 


‘that for nonaerated flow. This assumption was tested by making oe 
|e the same channel under nonaerated conditions and comparing the results. — 

| This comparison ‘showed that for ‘the rough channel used in the present experi- 
ments the resistance coefficient was essentially the same for both core ag oll oi 


_ With this agreement in resistance coefficients it was possible to ‘compute a 

the velocity and depth of nonaerated flow having the same water discharges Sd a 

and slopes as the aerated flows used in the experiments. . Figs. 16 and 17 show | “ 

the ratios of the aerated depths d, d and d 4 to the computed depth for an equi- ae 


valent nonaerated flow dn n aS well as the ‘ratio of the 1 mean 


Proc. Paper 1890, December, 1958, by Lorenz G. Straub Aivia 

Director, St. Anthony Falls Hydr. Lab., Prof. and Head, Dept. of Cw. 

Prof. of Hydromechanics, St. Anthony Falls Hydr. tab. Univ. of oe 
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- i the usual laws of open channel flow. Now if the channel is hydraulically = aaah 
rough so that th iction coefficient is independent of the mixture Revnolds 
* 
— 
* 
ue 
— 
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mean water velocity was expressed as 


in which V is the local velocity and C is the air concentration at elevation y ¥ 
velocity for an open n channel flow for the same and 
- slope was computed in the usual way in accordance with Eqs. (17) and (18) — 
a just as though no air were entrained. — _ ‘The ratio of this velocity ' Vi m and t the Ps 
3 _ mean velocity for the equivalent a air entrained flow was determined and plotted 
An Fig which is reproduced wi e new points added to it as Fig 
ee, Fig. 17, which duced with th ts added to it Fi od 
_— 2a agreement C of velocity ratio based ‘upon velocity measurements a and } 

based upon air ‘concentration measurements is ‘remarkably good. 
ae _ The computation of the normal nonaerated elements of open channel flow q 

ina particular application 1 forms a base from v which the elements pertaining x is 
te aerated flow depart as a consequence of the entrained air. ‘The data and ane ; 
curves given in Fig. 12, 16, and 17 form the basis for estimating the depths 
i a steep channel of equal roughness on which air en- | 


The value of the Manning n for | the » experimental PORES 


— depth d asd able to re wi to the 
discharge ibution curv ity in the a yield 
ation distr f the veloc flow i 
- tr surement 0 way a 
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; varies from ; about 0.012 to 0. 016, depending upon the depth of a ‘so at eet 
channel roughness corresponds approximately to a surface roughness 
such as may be found in practice. In a particular instance of flow on a steep Re 
channel the water discharge and slope are known from which S/q1/5, where S_ 


to the curve in Fig. 12 will supply ; an estimate of the expected mean air Sgn 
concentration. Entering Figs. 16 and 17 with this mean air concentration a 

_ yield the various depths and the mean velocity of the stabilized air entrained ae 
flow in terms of the equivalent depth me mean ehtestend — in the usual — 
way as and flow were no’ 
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a Federal Communication Commission allocation for the operation of 10-cm 
_ radar for weather purposes. Since 10-cm radar equipment is used by the bcs 5 
Federal Aviation Agency and the Air Defense Command for aircraft tracking, 
8 is reported that the F.C.C. plans to limit the use of weather radar to the 
The discussor questioned the use of rainfall data as a 
echo patterns. It should be pointed out that in this 


ye of ‘echo pattern during the hour. " ‘Therefore, it was felt that the use of S . ; 


the echo pattern at the midway point between the hour would be representative 


ample, ‘if Pr is 4.62 decibels higher than the equivalent power from the rain 
gage, the error is 50 per cent in rainfall amount regardless of the time, Se in 
whether it is ten minutes or ten hours. On the other hand, the integrated value 
of the difference in powers in decibels with respect to time will vary for at 7 
‘different storm lengths. It would also appear that a more natural base (cutoff i 
value) would be the minimum detectable receiver power. The writer sees no 


Treason to believe that this base should in any way be related to — 


. Paper 1901, January, “1959, by Glen E. Stout. cag 
3 “State Water Div., ‘Urbana, 


e author greatly appreciates the careful reviewand_— 
q concerning the subject. Since publication of the 
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_ WILLIAM H. SAMMONS,1 M. ASCE. ~The research reported by Mr. Benson 


peta of an optimum main channel- ‘slope factor for us use fe in flood- frequency 
analyses in New England. This generalized approach - may be expected to yield 
: improved channel-slope factors and is to be recommended when the nature of 
problem does not justify detailed study of channel- slope ey 
4 The paper is extremely stimulating on the history, | philosophy and the - 
_ situation of the treatment of hydrologic factors by generalized statistical © 
; models. . This presentation is timely and d valuable t because hydrology as ner 
a science or a technology has come of age and the statistical methods have 
been recognized (demonstrated) as an essential “tool” for the scientist and 
_ Owing to a general negligence. of the problem, the need for an lathdiidae a 
f basic- data program in hydrology should be recognized and must be made a 
— to all concerned directly or otherwise affected. The author has given Tas 
us some insight into a number of the hydrological ‘parameters entering into” = 
_ the prediction of peak discharges on a generalized basis. His short presen- is 


that all of the complex mathematical expressions concerned with the behavior 
‘of the hydrological cycle can be automatically and simultaneously solved i 

simply by pushing a few buttons - Although something of this sort may eventu- 
ota come to pass, the conversion of the complexities of the hydrological — 


| great m many years. Much anes has yet to be developed. In the absence of i 

tested academic theory, it is customary to set up “models” expressing the _ 

general characteristics of hydrological situations. A model usually has 
basis in theory that has been “condensed” to certain dominant factors or indi- — 
ces that can be evaluated by means of field observations of representative 2. re 

“length of record, recent maps and surveys, etc. Ref 

The author states that “ “The main~ -channel has been found next in 


ept. 
Beltsville, Md. nial Unit, U. 8. Dep 
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"England that ‘ “other” hydraulic factors, yet unknown or untried, could also 
qualify. ‘Itis necessary to examine basic assumptions(1, 2,3,4 5) of logarithmic 
ae regression models as employed by the author. The logarithms of the original 
i Ee variates were assumed to be transformed to a log- ~normal probability distri- ‘| 
5 i; bution. This may be true to a certain extent p provided the ‘original variates — 
all have the same or similar probability distribution (characteristics). The : 
atter depends on many items, some are known, others unknown. c An example 
— to illustrate a similar situation but converse to the author’s New England find- 
—_ is presented. A log-regression model between median seasonal yield, _ 
(1,000 ac-ft. units) as the dependent variate, median seasonal watershed 
precipitation, Po, (inches), and drainage area, A, (square miles) as the inde- | 
pendent variates should illustrate the pertinent facts and modifications. | Con- | 
sider the variates tabulated as shown in Table 1- ~a, “Results o of ‘First Approxi- 
Ee ‘mation’ ’ Multiple Regression with Yield”, with their respective statistics © 
(first approximation). A, graphical examination of the cumulative probability 
_ distributions of the original variates indicated the need to rectify the variates 
_ for normalizing. See Table 2, “Statistics for Original and Transformed Vari- 
ates with Yield”. From the statistics tabulated in Table 2, rectification = 


6) to normalize the variates were by the 


: 


ent of variation of the original variate 


corn cv is the coefficient of variation of the transformed variate 


ist the arithmetic mean of the original variate bid 


ae An empirical rectification of Vg, A and Py would be (V2 + 124), (A tum 

and (Pg + 68) where = 124 (1,000 ac. ft.); =1172 (sq. mi.); and 
= 68 (inches). — The rectified variates gave reasonable cumulative atta 

bility distributions fe» re-plotted on log- -normal probability paper. he! 
The new statistics (second approximation) are tabulated in Table 3, “Re- aa 

sults of ‘Second Approximation’ Multiple Regression with Yield”. A compari- 
4 son of the statistics tabulated in Tables 1-a and 3 indicates that a curvilinear — 
logarithmic model would best represent the data. (The author’s Fig. 3 shows 

a similar trend because of the curvilinear regression constants of a etc.). ose 

_ For all practical applications the following simplified equation would suffice: 


ere A is the in square miles 
Vo is the median yield in 1,000 ac acre-feet ‘unite 
_ Tables 1- “a and 3 also indicate that precipitation was more important than 
the drainage area, which | is contrary to the author’s findings in New England. 
_ Table 1-b, “Results of Multiple Regression with Peak Discharge ” confirms the 
=. statement which was indicated by Tables 1-a and 3. In the latter case, 


‘median 2 annual | peak discharge, , Q2 is related to A and P2 in the same region 
(an additional gaging station is included with the peak discharge relationship. 
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Table 1- b. Results of Multiple Regression with Peak Discharge 
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Coefficients & Estimate "Student"- 
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Table 4, Statistics 

Constant: “Reg. Coef. 
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_*ZK will be referred to as Z-values for short. ae if 
a-values were "adjusted" "(smoothed by eye) 


: - The b- values for T = 50- and 300-year were omitted — 
non-representative in the ‘derivation of the above relation 


‘ship. 


for New England ‘and Relation Between Selected Recurrence 


where etendaré error of was 
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i Table 3. Results of "Second Approximation" Multiple Regression with ris | off 
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region of the western United States where such results are to be meee. > g 
Because of interactions(7,8,9) between variates, it is very important that | 
all possible combinations of the variates be investigated to determine the a 
- efficient as well as the most logical order of introduction into the regression — 
a _ The author’s statement Coupe provides the best correlation with flood magni. a 
tudes” is at first misleading until Figs. 1 and 2 are examined and thoroughly : 
"understood. _ A formal definition of the standard error would be of assistance | 
to the readers without a statistical background. The variation of the ‘standard 
error with the s slope factor is the “key” to understanding the figures. 
_ The term “generalized” versus “regionalized” has been used in the litera- 
| ture with interchangeable meanings. . The wi writer riter questions which is correct a 
| The author refers to the New England study as a “first phase” and continues 
to state that runoff, precipitation, topographic quadeonaie maps and historical 
flood data are available in this area. The writer wishes to emphasize that 
"studies of this type must precede an ancienle and be readily available. The 
meteorological data should be continuous from region to region. It must be of 
sufficient ¢ detail to indicate sub-area details, otherwise the variability which — : 
must be explained would be lacking on the maps. An areal integration of the 
normal annual precipitation (or runoff) maps should show a normal (or log- a - 
normal) probability distribution of the specific variates, otherwise the 
aly variate should be rectified prior to the introduction into the regression 
model. In order to realize the above, the Hydrologic Investigations Atlas, Om 
| ae as HA-7, should be extended to cover the continental United States (and 
parts of Canada) | future studies may be realized. 
_ Topographic quadrangle maps of a uniform standard a are lacking. © Some cand 
states have very recent coverage. . The adjacent states (within the hydro ological — ‘f 
region) may be in the “horse and buggy” state of 1 mapping. Organization in r 
. | the government agencies and local governments is inadequate to effectively oe 
solve this problem. Topographic parameters relate “only” that which 
been recorded in our present maps by s states and not by hydrological eagiees. 


region whereas the same index will fail to chow significance in the adjacent “a 3 
region due to lack of uniform standards for collection of field data. A peibers 
| program to fulfill the future needs of hydrological data is desirable oF Wee itt a: 
Historical flood data is indispensable in any regional analysis. ‘The data 
‘indicates and may give ve guidance not readily apparent in the raw data. Caution =) 5 
should be exercised when decisions are made based solely on historical ldata, ak 
‘ioe experienced analyst, like a good de detective, may wield | the historical l facts — 


e: Each variate(5) should be be presented as a histogram, frequency polygon or 

cumulative probability distribution. The original (observed) and 

| (rectified, etc.) data may be visually compared (' (the cumulative probability — a: 
distribution has many advantages o over the former) . The author deals with = } 
| ranges, etc. which are difficult to visualize. 2 
'a The statement “. . . the use of mathematically fitted pa involving prior 


as to distribution might obscure the v relationships that were 
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possible to the crtginat data at nile station on would b be most desirable” is 
a weakest link in the regional type investigations, in the opinion of the writer oF & 
fe A priori distribution selected in a haphazard manner or by a policy decision | 


is very undesirable. Mathematically f fitted curves to the “ “proper” ” cumulative i 


probability distribution will eliminate much of the bias and errors of intui- cs 2 
a 4 _ tion involved in graphical analysis. Proper rectification (end transformations) | 


of gaging station observed data will eliminate the bias in graphically drawn | 
ra curves. Plotting of the ordered observed data is of minor importance provid- 
_ ed the analyst is consistent and does not incorporate fallacious interpretations 


the results. selected for | the ‘regional ‘Study may dictate a priori 
of the variat 


(as for a in Fig. 3), thus the linear trend of sae regression coefficients b and 
i c in Fig. 3 (log-extreme value probability paper). The graphically construct- 
Ra 5 ed curves and the distribution of the variates at selected recurrence intervals 


i creates the concave t upward bias of the constant term a on Fig. 3. ‘silbigedd om 


ae Plotting of the constant term, a, and regression coefficient, b, for QvsA 


is leads to linear relationships. _ (This point will be demonstrated in an example : 
in the latter part: of this discussion) 


mathematically fitted curve dealing with the basic statistic 
- the distributions would involve only one, two, or possibly three ssiitehe. cal 
_ correlations of the variates. B. The latter approach(10, 11,12,13,14) has been _ 


a _ and agencies. Limitations of the basic data (observed or original) - 


ae followed by the writer and ‘many other statistical hydrologists in various uni- j ‘g 


tend to support the latter approach rather than the author’ 's refinement. _Sys- | 
tematic curvature of the cumulative probability distributions on “Type A” 
} _ probability paper suggests the use of “Type B” probability paper, thus linear- 


izing the data with the proper transformation. 


Cry 
_ The author presented eight parameters which were studied. The writer _ 


would like to point out that in all cases where the “average” was used a nil 


_ substitute of the “median” would prove as good or less variable and more 


a 
satisfactory in a regional analysis. . The of the average rather than the 
a, —aae can, and usually does, introduce uncontrolled variation which is very 


difficult to explain or rectify. In summary, the cumulative probability distri- | _ 
bution of each prospective parameter should be known (not assumed). Weight- _— 
‘“ ed observed data should be supported or justified before the use of indiscrimi- 

nant weighting coefficients which may only add “insult to injury” in our re- 

a The use of topographic maps to obtain elevations and distances for — 7 
; strani profiles depends on the scale ratio. Gaging stations are known points 
(elevation) but considerable error can exist in the profile elevations obtained | 
between gages in rugged topography. a Small watersheds 1 require tolerance 

- standards which can not be defined by the present United States Geological — 
Survey maps (scale - att = 2000 ft. = ay 1/2 minute quadrangles) and must be } 

_ supplemented by er engineering field surveys. Did the author make a compari- 

er son for profiles obtained from different map scale ratios versus a field - 
survey?” If so, what are his recommendations? What difference, if any, 
be noted or the determination of the best slope - factor by 


ies 
HY 
em 
col 
Fic 
8 
COE 
3 int 
my 
to « 
okt | val 
pro 
abo 
ee. 
a fitt 
= 
ae 
h 
a 
To 
) 
a or t 
4 
be 
4 e 
T 
tl 
ul 
an 
‘ 
fc 
 ilations for a specific plotted point on Fig. 1 would clarify the technique Felat 


employed (the: author that “... the 10° slope is far simpler to 
ce The author’s Fig. 3, “Variation of Regression Coefficients with Size of euler om 

” is very informative although elusive to the casual reader. The writer 
contacted the author(15) and obtained the a and b coefficients for recurrence 
intervals, T = 1.2, 2.0, 2 1/3, 5, 10, 25, 50, 100, 150, 200 and 300-year events 
to demonstrate one of the many useful relationships that may be developed. — 
Additional information on the hydrology of the region is the result. | BPs, 
- “Statistics for New England and Relation Between Selected Recurrence Inter- i 
| vals and Drainage Area” tabulates basic regional information. = = = =~ 
‘The a-values plotted as a smooth curve with a concave downward trendon = 
probability paper similar to. the author’ s Fig. 3._ ‘The 50-year, a- -value plotted ool 
_ above the curve and the 300-year below. " Entries in Table 4, Column 2 are 


the b-values (unadjusted) were plotted on similar probability paper ‘and on ua 
7 log-log paper . The 50- and 300-year events were also non-representative ‘aaa 
and were omitted from the computation of the trend line. was 
“fitted with the following best fit equation: 


= 0.889532 


0088 : is s the standard error of e of ‘the 


bs ‘The only reason for selecting the log-log plot over the extreme value log- 
probability plot was to simplify the computations. An the latte er 

would take the following form: 


where log 1/a is the 

7 


: as requires a table of y-values for selected recurrence intervals. ar 


The  y-values are not in wide use although tables(14,16) have been published. — 
re- | The computed b-values from Eq. (3) are tabulated in Table 4, Column er <: da 
ae _ The assumption was made that since a log-regression model was aaapreyes 
lotting by the author, the log- ~normal between recurrence intervals 


bility distribution. From Column 5 of Table 4, Table 5 was aceapelal aah 
tabulates tRo- -values for T and A. See Table 5, “Relation Between —— 
Recurrence Intervals and Drainage Area.” A plot of thistable on log-log _—~— 
paper gives a graph of expected variation hetween the magnitude of selected ‘fh 
intervals in n relation to the median (2-year event) and 
Table 6, “Relation Between T, Cy and A, New England” was 
| the basic(1, 2,3, 4 17,18 19) log-n -normal | distribution | 
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table of ZK-values is in References 3 and 


is the coefficient of variation on of the ‘original variate dis- 


cussed in References 1, 2, 3, 4 and 


Drainage Area in New England (Peak Discharge) 
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_ The following 5 cumulative probability distribution iia 5) of Cy v a 
Tare given in lieu of a family of straight lines on log-normal cednihier ot 
| paper for A = 1, 10, 100, 1000 and 10,000 be ph -mile drainage areas: Jie 

+0. 1901 K for K for = 0.26590 bee 


Gy = 006880 0.0667 K for = 0. 


= = 0. + 0. 0272 K for = 
Baht 
= 0.54K25 + 0. K for C= 
‘Table 6. Relation Between Recurrence ‘Coefficient of Variation 
and Drainage Area in New England, Peak Momarge 


Baste Date from M. As 


Ke. 
Average 1.1465 0.8976 0.6 0.5074 


i T is the recurrence interva 


= 
cy 0.8976 + 0.1188 K for ¢ = 0.1326 


t 
(5) 
(6) 
— 4 
d 
= 
i — 

— 
al 
"ome 
— 
the areinage aren in square miles. 


ae the K-values for specific levels of probability must st be ‘selected from % a ; 
each published in References 2, 3, 18 and 19 for the Cy-value given following 


years of record” has been employed in a log- 
ra gaging station. 


“Student's at the 10% 1 level of ‘significance, ¥ with 


aa vent Eq. (12) and Table § 5 the following tabulation was vas made: 
Drainage Area in Square Miles 


1000 


_:100-year recurrence interval peak discharge, and if the ‘ratios, 100R2, should 
a ‘used to evaluate the expected minimum years of record necessary by com- 


ee puting Ym from Eq. (12), it then could be concluded that 57 weare of record (or 


i = area whereas only 13 ye: years s will be required for a 10, 000 square i 
mile watershed. How many watersheds of 1 square mile have 57 years of _ 
record in the New England region? Needless to say, drainage area alone is 


= not sufficient in New England to estimate peak | discharge. The addition of the 
eae author’s slope parameter (and others) should decrease the minimum required 
_*years: of record, Ym = = 57, to some smaller more = eee | for al- aaa 


watershed. 


Variates for the Logarithmic- 1958. (Mimeo.) 


Chow, Ven T. , Discussion of of “A Log og-Normal Frequency 
Analysis”, -TAGU, Vol. 38, pp. 962-963, , 1957. 


ae Kleen, M. H. and |W. H. Sammons, “ Agricultural Watershed Influence on 
foe: Runoff Potential” , Manuscript presented | at the thirty- ninth annual meeting 
c of the American. Geophysical Union, May 5- 8, 1958 at Washington, D. C. ae 
May 1 1958 (Mimeographed). ke 


drainage area is considered in the New England region to estimate the 2- i ta 
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@é ji. 1. Aitchison, J. and J.A.C. Brown, “The Log-Normal Distribution’, 
— _ Cambridge at the University Press, 1957. 
Brakensiek, D. L., “Fitting a Generalized Log-Normal Distribution to 
— 


. Goodrich, R. Dra Line Plotting of Skew — Data”, 
Valavanis, ‘Stefen, “Must the Diffusion Index Lead? ” The 
Vol. 11, No. 4, pp. 12-16, October 1957. 


8. Maher, John E., “Correlation Among Time Series and Moore’s Statistical — 


Indicators” The American Statistician, Vol. No. 4, pp. October 
9. . Moore, Comment on “Must the Diffusion Index Lead?”, 
‘American ‘Statistician, Vol. No. pp. 16-17, 1957. 
aa Schwarz, H. E., “Determination of Flood Frequencies i in a Major ‘Drainage a 
nae Basin”, Washington District, t, Corps of pe of Engineers, April 1957. (Mimeo- fg 
Corps of U. S. Army, “Stream Flow Volume Duration Frequen- 
, U. 8. Army Engineer District, Washington, June 
. L., “Statistical Evaluation of Runoff Volume Frequencies”, 7 
Sammons, W. H., 7, “An. Evaluation of Procedures for Development of Short 
Rainfall Depth-Return Period Relations from Records of 


recording Rain Gages”, Thesis, Grad. School, Univ. of 
. West, E. M., and W. H. Sammons, “A Study of Runoff eee sia Drainage 
Areas and the Opening in Attendant Drainage Structures”, Report No. 2, 
Highway Materials Research — of 


15. . Personal ‘Communication with M. J ‘Benson, ‘April ‘17, (1959. 
“Probability Tables for the Analysis | of Extreme- Value Data”, 

Weiss, L. L., “A Nomogram for Log-Normal Frequency Analysis® 
TAGU, Vol. 38, No. 1, pp. 33-37, February 1957. 

Chow, Ven Te, “The Log- Probability Law and Its Engineering Appli- 


cations”, Proceedings 80, Separate No. 536, November 1954. 
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ee Discussion by William C. Peterson — 
| WILLIAM C. PETERSON. 1—In computations of water- surface profiles by = 
“the step method, it is often necessary to adjust a velocity head computed from 
the average velocity in a subdivided section to a value more representative of 
the true velocity head. _ Adjustment is usually effected by the use of an ap- - 
‘proximation of the correction coefficient 


= 


ties by stipulating (as do the authors) that in a matin 4 section 


= 


then 
s With « o< expressed in terms of section properties, the computation of oc , 

_ can be incorporated in the authors’ Table 1 and eliminated from involvement 7 7 a | 
in the trial-and-error step procedures. _A revised format of Table 1—from a — 
_ which curves of A, Kg, , and oc o< versus elevation can be prepared— is shown in — 


iS In the authors’ Method A the arithmetic mean of the friction gradients at 
_ the ends of a two-section reach is used to obtain the friction head, hy. This ; 
| = & use of the arithmetic mean overlooks an interesting re relationship, ‘inherent in 


: a. Proc. Pape: Proc. Paper 1997, April, 1959, by Joe M. Lara and K. B. httiennie.. 
Hydr. Ragr., ‘U.S. - Geological Survey, Los Angeles, 
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that becomes: evident when the is rewritten 


: 


“That is, Q is a function of the qoematvte: mean of the end- -section acl 
a and the geometric mean of the square roots of the end-section friction 
Expressing the friction gradient, Sf., in terms of a “Project Q” and the 
: conveyance, Kg., leads to a step method of determining water-surface pro- 
files that requires only the 12 columns of tabular data shown in Table 5. ~ te In 
this method trial values of water-surface elevations are successively as | 
sumed and the unique value is confirmed when the “Computed Q” equals the 
_ “Project Q.” Columns 3 to 6 of Table 5 provide data to determine hg i in 
column 7. The square root of the friction gradient, St.., is obtained in column 
9 and used with the conveyance, Kq,, to compute Qincolumn12. 
Bs _ The formula for hg in column 7 follows from the authors’ Fig. 1: 


— 
| 
1 q 


DISCUSSIO} SION 


Assumed 

on water 
surface 


Section 2} 


1. By letting (1- ‘M) = E-= = an energy-conversion factor, the 


is 1 to contraction expanding reaches provided complete 
 (E= 1.00) conversion of potential head is assumed at all times for contracting : 

. reaches. ‘This assumption—equivalent to assigning a value of zero to M for 
contracting reaches— —appears reasonable when one considers that the total — ai 
angle describing a contracting reach is generally much smaller than the maxi- 

_ mum angle of 12.5° in converging transition structures for which a value of _ 
0.1 is usually recommended(2) for M. For expanding reaches (downstream hy 
senainty less than upstream h y) the value assigned to E may range from 0.0 — 

As “hydraulic characteristics” were made available only for Section 1 in — 
‘Table 2 (Method A) of the authors’ paper, the sample computations in Table 5 — 
are limited to those for a trial water-surface elevation at Section 2. Needless — a 
to say, the ; assumption of no separation losses in this contracting reach and 


the use of the weighted conveyance, Kq,, will result in a slightly different _ 
"Method 


Computation of the conditions at crossings take ad- 

7 _ Vantage of the wealth of experimental data obtained by Kindsvater, (3) ae a 
* - Carter, (3 4) and Tracy,(4) in their studies of open- -channel constrictions. The 
a. “bridge- -crossing” computation would be handled best, perhaps, as a separate a 

_ investigation and the results integrated with those of the step procedure. — Siteene. 

. ‘The authors’ ’ Method B makes provision for a condition where the length of 

- the overbank flow differs materially from the length of the channelized part _ 
of the waterway. Undoubtedly, field verification of Method B exists — 

warrants its use in preference | to Method A. writer welcome 


gem utation of Peak Dischar; e, administrative Department of ofthe 
— U. S. Geological Survey, 1948. 
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Hunter | Rouse, John Wiley & Inc., New 
“Tranquil Flow Open Constrictions,” by CarlE. | 
_ Kindsvater and Rolland W. Carter, Transactions, ASCE, Vol. 120, 1955. ie ail 
4. “Backwater Effects of Open- -Channel Constrictions,” by So : 


Transactions, , ASCE, Vol. 120, 1955. 
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by Peter Ackers 


ae PETER ACKERS. 1—The number of papers published in this journal in 
: recent years on the resistance effects of various types and patterns of rough- a 
ness shows the great interest of this subject at the present time. The authors’ 
Work on large- -scale triangular flumes contains very useful information, which 


has particular significance because the possible effects of the variation of the 7 


shape of the flow section with depth are eliminated with an open channel _ ff 
this form. rri¢ 
a roughness experiments on n salt- -glazed : sewer pipes with various eccentricities | | 
imposed at their spigot-socket joints; and, contrary to recently- -expressed 
3 - opinion, (1) he found that the Colebrook- -White(2) equation (on which the Moody ; 
diagram is based) agreed satisfactorily with the experimental results over 
¢ the full practical range of operation. (3 
soma Joint spacings of 2 ft and 3 ft were used, and eccentricities varying up to ad 
_ 0.4 in. were imposed at the joints. The data was used to evaluate €, the 
equivalent sand roughness in the equation used, and an attempt was ‘made to 4 
correlate this with the joint pitch and eccentricity. | So as to retain dimension- _ 
al equivalence, the experimental value of € for each roughness condition was 
al compared with h2/) (h being the mean height of the lips \ of the joints and A ——_ 


‘ their pitch), and an approximately linear relationship was found to exist be A 
wom these variables. With this in mind as a possible basis of correlation, a i 


‘Posey, which may be of interest . The authors had reported one quite unex- — ral 

_ pected result; that without battens their triangular channel appeared to be oe. i 
“smoother ' than smooth’, and i they had to modify the usual smooth-pipe equation — 

_ for this reason. In consequence, the writer also had to modify the Colebrook- 


ter aleo e Colet 


effect, arising from the non- 1-uniformity of boundary ‘shear around the perime- a 
_ ter of an open-channel, and the author’s opinion on this would be welcome. oa a 
ae Confining attention firstly to tranquil flow (Froude number <0.7), mean _ * x, 

values of € were calculated for each batten spacing, and the ‘suitability of the " 


Prin. Scientific Officer, ‘Research Station, ‘Berkshire, 


NTS IN A TRIANGULAR CHANNEL? 
| 
| 
q 


0.0539 ft : 


4 channel, neither more nor than the authors obtained with the 
7 _ Manning equation and their modified Prandtl formula. With battens at 24 inch | 
orm, the Colebrook-White equation is a good fit to the experimental data, a = 
but it deviates progressively as the isolation of the Toughness | elements de- ‘ 
creases. The implication is that the shape e of the transition between ‘smooth- 
a , turbulent and rough-turbulent flow varies with the batten spacing, but as open . 
& channel data does not yield lines of constant relative roughness on the usual 
- friction factor versus Reynolds ‘number diagram, this can not easily be illus- 
33 trated. However, the transition would probably change progressively from a a w 
descending curve at high values of pitch to a rising curve, such as occurs 
with corrugated pipes, 4) y with closely spaced battens. ." It is interesting tc to — f 
serve that, in spite of the short-comings of the Colebrook-White equation 
which | are apparent with this type of roughness a at close spacings, the 
ship between the € and values | given above is roughly linear. However, 
it would be unwise to generalise on this basis, as the available data covers -_ 


batten n height on only. 4 Do the authors plan to extend their ir experiments to in- 
writer has plotted the equivalent roughness against Froude number 
‘ig. 1, confirming the author’s report that resistance is | higher under 
critical conditions than for tranquil flow. . Although there isa great dealof _ 
Am - scatter, the plot of the data shows a trend for the apparent roughness to be ~~ 
q greatest at Froude numbers of about unity, when surface effects due to channel 


=) 


M. Morris, “Design methods for flow in conduits” , Proc. ASC 
HY 7, July 1959, p. 43, a 


* « F. Colebrook, “Turbulent flow in pipes, with particular reference to the = 


“3 transition between the pipe laws”, J. civ. Engrs., 


Anthony Falls Hydraulic > Laboratory, Technical paper No. 
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RESISTANCE PROPERTIES OF SEDIMENT-Li LADEN 


- Discussion by Emmett M. Laursen 
‘EMMETT M. “LAURSEN, M. ASCE. —In this paper the authors very | 
demonstrate the truth of their statement that . . when a stream has a 
able bed, and sediment is being transported, the problem of determining the 
_ resistance is much more complicated than in the simple case of clear water 
_ flowing in a channel with fixed walls.” On the basis of their observations and > 
3 ‘measurements they conclude that the interaction of the sediment and the water 
as it effects the resistance to flow can be considered in two distinct ways. _ 
k The larger effect is due to the bed configuration (whether ripples, dunes, or 
- flat)— the conditions of flow, and the properties of the fluid and the sediment : 
Ee: in some still unknown manner creating a certain bed configuration, and the © : 


roughness of this bed configuration affecting the resistance to flow. The ob- 
served qualitative relationship between roughness and resistance is entirely 
in accord with experience with fixed oe and the magnitude of the orn is = 


The smaller effect the authors is. a decrease in 

resistance to flow with an increase in sediment transport—all other conditions - 
. - being presumably unchanged. This effect they ascribe to the damping of the * 

turbulence by the presence of ‘suspended sediment. Certainly, the presence > 

_ of the suspended sediment must effect the flow and it may well have the | effect oh 

on resistance suggested by the authors. | 
& However, the writer would like to suggest that this is not yet proved beyond ox 
= reasonable doubt since the sediment load could conceivably effect the re- “ie : 


in another way. Considering a flat bed for simplicity, there is 


-— is added to the flow over the fixed granular bed (returning to the con- a 4 
3 ditions of the unfixed bed). In the case of the flow of clear water over the . Qa vs 


fixed granular bed t there is a force transmitted to the peg ; through the pressure . 


added, there will be bed load movement as as load movement. _ 
‘The resisting force of the fixed grains are now transmitted to the flow at least least . 
partly in a different manner. particles rolling and sliding along the 
boundary as bed load will be moving at a velocity less than that of the fluid 
them in order to receive a propelling force from the fluid. 


q 


= 
| 
| 
3 
3 
Associate Prof of Ci by Vito A. Vanoni and George N. Nomicos. 
of Civ. Eng., Mich. State Univ., East Lansi 


the fixed boundary should be somewhat different anes on whether or not | 
the water is transporting : a sediment load. 
ie Whether or not the rolling and sliding particles of the bed load ‘movement _ 
_ could account for the observed difference in ‘resistance could be checked ag 


| 


— 
_ vember 1080 a 

q 
Se 

| 

“gt 
| 

i 


KET REQUIRED TO PREVENT WAVE EROSION? 


"Discussion by P. Bruun 


BRUUN, 1 F. ASCE.—This and to those who 
are trying to introduce better construction practice within this particular 
z field. Meanwhile, it includes two important assumptions which may not a 
ways be easy to keep up with in practice. One is that skilled labor is avail- 7 
&& able for such construction work which requires great care. . The other is the 


problem of availability of | proper material for a blanket (or rubble mound) . 


In regard t to the first question, it is probably true all over the world -aughedl it 
Btn _ One essential reason for this difficulty is the replacement of hand-— 4 
work by work done with the a aid of mechanized equipment. _ This makes possi- 


; the requirements of accuracy because connection between the brain and hand» 


& man is closer than the connection between the brain of man and the oper- — 
ation of a mass-producing machine which can never be as flexible as the hand- . 
tool. The question of skilled labor has sometimes become so acute thata 


: different design had to be worked out. Examples of this are the replacements _ 


"with more rugged and less detailed designs. Replacements range from use of 
; unsorted quarry waste material in thick layers under the cover stones (usual- 


ly making the design less stable thus requiring more repairs) to monolithic 
ss In other cases lack of proper filter material—within economic limits—_ ~ 
prohibited the use of gravel filters and this circumstance changed the design. ZZ 
_ Florida and the low countries of Europe present many examples of designs of 
- jetties, groins and seawalls where the material available was the determining — 
x factor in the final design. - Attempts have been made to replace gravel filter — 
layers and rubble mounds with devices while still preserving some of the good . 
qualities of the normal graded filter layer of traditional a ee 
_ Willow- -mattresses were used for decades for bottom protection against | 
scour which means as filter layers for special purposes. Nowadays these Bs - 
= been replaced, among other measures, by -asphalt- matresses such as 
coastal pro-- 
tection and harbor works). . The most modern development is represented by : 
4 nylon sheets which are now being tested (in common use) where the bottom is 
to be protected against scour caused by : strong current activity e. .g- where pate 
caissons for dams, sluices, weirs, and jetties are to be built. Dutch 
Proc: Pape Proc. 2021, 1959, Enos J. Carlson. = 
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_—-Bublication on the “Delta Project”, July 1959, No. 8 ‘shows how 3 m 

thick nylon sheets can be used. The main difficulty in the application of those 


sheets still is the technique of bringing the sheets safely to the bottom of the 


ee nylon sheets are entirely impermeable for sand and impermeability 
% 4 desirable for the work the sheets are supposed to do. . Meanwhile the true at 
- definition of a filter is still a layer which, while letting water pass through, a 
holds the material back. It is undoubtedly possible to use permeable she sheets - 
the material to be protected is not too fine. 
_ About two years ago a member of the staff of the Coastal Engineering Labo- % 
"ratory of the University of Florida, J. J. Leendertse, got the idea that Fiber- 
_ glas mat material of insulation type might be useful as replacement for 4 
filter layer in a Florida seawall. The laboratory contracted proper Fiberglas 
manufacturers who furnished mat material free of charge for field experi- ts 
= ments with Fiberglas as underlayer for 190 lb. interlocking concrete blocks 
in a revetment resting directly on a sand slope. _A plastic cloth manufacturing» 


_ company later asked the laboratory if their plastic woven material would be 2 


2»: 


useful in . coastal structures. ‘The laboratory found some of its material suit 
Af = and other field experiments with this material were started at a lower 
_ elevation (around sea level) than the mat material. Later the proper interests 
- entered into a contract with the laboratory on experiments with | the application 
- of mat as well as cloth material (Fiberglas, plastic mat and woven material) © 
: as filter layers in coastal structures including rubble mound jetties, , Seawalls, 
‘revetments fc for coastal protection, harbor 1 works, vertical bulkheads, 
 sheet- -pilings. The result of field and laboratory tests which are now being 
_ published by the Engineering and Industrial Experiment Station of the Uni- - 
Bibenwe sy Asi: of Florida demonstrated definitely that Fiberglas mat material with 
some improvements now being made is well fitted for application as filter 
P 7 _ layer above sea level in several coastal structures, particularly revetments, 
oa bulkheads and seawalls. The mat material was, thanks to its loose structure, 


; a applicable » below sea | level where the cloth material proved | successful 4 


particularly in rubble mounds, protective aprons and other bottom protection 


: -* The same cloth material seems to be useful for tightening leakages” 
_ in rubble mound jetties (replacing a asphalt- ~grouting) and it may prove see 


ag 


‘tificial no ment may be increased. = 


Re _ ‘The » intention of these remarks are in : agreement with Carison’s s desire to 
— stress the importance of well built filter layers of traditional type as well as wk 
inventions. _ Carlson’ this to attention much credit. 
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Discussion by Robert O. Thomas. 
ROBERT O. ‘THOMAS, 1 F. ASCE. —Ever- -increasing demands for additional » 
. water supply development ; are resulting in a deepened interest and appreciation 
of | the utility ¢ of ground water storage for conserving, regulating, and distribut- — os 
4 ing an appreciable portion of the quantities of water required for maintenance * 
of our expanding standard of living. _ The following discussion stems from a a 
a more or less close association with | much of the thought and study for the de- a 
velopment of ground water resources in California. 
Conclusions as presented in the paper concisely point out that the existence 7 
_ or occurrence of water in the interstices of the materials composing the wal 
crust of the earth is not a a separable, independent physical phenomenon but is eo 4 
inextricably linked with the preceding and succeeding phases of the hydrologic — ; 
cycle. In passing through that cycle fr from precipitation, through infiltration wi 4 
and surface runoff, to capture and use, either by nature or by works of man, -<_ 
and finally to evaporation and return to the atmosphere, to fall again as pre- 
- cipitation, the water droplets many times pass through the state of temporary _ 
— in underground storage. Such a broad concept is —— ; 
_correct when made without regard to time, whether such be counted in mplautes a 
or in eons. Consequently, all subsurface water, even that which hasbeen 
trapped un underground for untold ages, isa part of the common supply and 
= to the natural laws covering hydrologic phenomena. It follows there- _— + 


from that man cannot increase the common si supply, although by the application . 
_of scientific principles to the conservation and use of the available water, he 
can so manage the available supplies as to increase the benefits to be gained = q 


through utilization of this renewed, and most important, 

- Utilization of ground water reservoirs for the maximum benefit to the area area 
ay involved is dependent upon adequate geologic and hydrologic data , from ce d i 
_ the results to be expected as a consequence of development can be estimated. = 
_ Prerequisite to the full development of ground water storage is adequate = : 1 
:* technical information, including the basic data necessary for application of ‘ee a4 
_ proved hydrologic principles. Required data can be secured by a continuing» a 
program for measurement of the amount of water passing through each phase a 
of the hydrologic cycle. With adequate hydrologic data, reliable estimates of — 


the quantity of water that can be made available on a firm annual basis by - a " 


coordinated or conjunctive operation of surface and subsurface storage or, 


er alternatively, by operation of the ground water basin alone, can be formulated. Pe 


a. Proc. Paper 2056, June, 1959, by Joseph E. Upson. ~ 


Engr. State Dept. of Water Resources, Sacramento, Calif. 
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of physical situations water will generally 


re the other in a complex pattern of skills and knowledge. “ One such area of re- 
& search may be be termed the geological phase, in which the conditions to be en- © 


; ‘countered on and under the surface are determined. The primary requisite — P 

is the location of underground storage space, which ; may be found in pomemens: 
Z formations that contain inherent voids in sufficient volume to hold, transmit, 
a and release water in 1 the quantities ne necessary to satisfy the demand on the _ 


ste 


supply. Such voids may be solution channels, as in limestone formations; oh 
-brecciated or fault z zones as in structural unconformities; or interstices be- — 
val tween the grains composing the formation. The: latter type is ordinarily found 
“ in alluvial valley fills a and cones and is composed of gravelor sand,or ~ 


ble material. . After determination that usable storage ‘space for 
a operation exists, the porosity, transmissibility, chemical composition, lo- Bey 
cation of intake areas, and other basic physical characteristics pertaining to 
hydrologic and hydraulic operations are determined. Mile 
_ The second field of investigation of subsurface storage areas is the hydro- 
logic phase. This phase is concerned with determinations of water Supply 
- available for the project; the regimen of its occurrence; the quality « of the a 
7 g native . imported, and subsurface waters; the operation of surface facilities 
for storing, transporting, and percolating surplus supplies in permeable | 
areas; the recovery of water placed in subsurface storage; and innumerable 
additional factors bearing on the entire problem of 7 
the place of origin and place of final disposal 4 
a: | economic phase is the third broad field of investigation. ‘This phase is j 
i _ concerned with the determination of places and amounts of use. It involves ey 
the classification of land areas as to suitability for irrigated agriculture, Bis “| 5 
domestic, industrial, and other uses; determination of water requirements for 
the various uses; production, distribution, salvage, and | disposal of water 
ia % ‘It should be understood that the above grouping of sie of water ‘supply - 
investigation is made solely for convenience in illustration and bears no re- — 
lation to activities ‘in various fields of endeavor. . The engineer and the geolo- 
gist will function in close coordination throughout the investigation and will be 
assisted as required by chemists, biologists, physiologists, economists, , botan- © 
ists, attorneys, meteorologists, and other concerned specialists. 
utilization of the ground water reservoir requires continual 4 
measurement of items such as quantity in storage, recharge, and ns a i 
a and, in addition, the rate of movement of water through the reservoir. The _ 
_ techniques for these measurements are more complex and the 
? Ag information needed to translate th the basic data into volumes of water is far i 
oes imposing than the requirements for a surface reservoir. r. Delineation of 
a areal extent and thickness of the rock materials of a ground water reser- 


sd is primarily a geologic problem. The quantity « of water in those materials ls | 
is dependent upon the porosity. . Permeability | is a major factor in the rate of = 
movement of water through the reservoir and therefore in the yield of wells - ; 
or springs. The geologic and | hydrologic characteristics of aquifers, and de- 4 


— 


* lineation of the water | table or piezometric surface throughout the pape water 
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Thereafter, in or 


Increments to storage also affect the position of the water table in re- 
_ charge area. Therefore, coordinated analyses of records of infiltration, soil 
moisture, storage and discharge, stream seepage, and ground water storage | 
in the r recharge area are needed for a thorough determination of ground water 
recharge. Similarly, determinations of natural discharge from the ground 
_ water reservoir require adequate data as to evapo-transpiration, outflow in 

4 springs and seepage to streams, loss of storage | by subsurface outflow pand = 
changes in the volume of water in storage. 

x Many reservoir problems are chiefly problems of movement of — with- 
in an aquifer rather than of replenishment to the aquifer. _ Developed aquifers : 

: serve both as reservoirs to hold water in temporary storage, and as trans- — 
_mission channels to carry water to ) withdrawal areas from the areas of re- 
charge. _ The perennial yield of a well or group of wells is determined by the 
quantity of water that can move through the aquifer from intake area. If the 

_ transmissibility of the aquifer is inadequate, the water levels in) wells Is will de- 
_ cline, whether or not the aquifer as a whole is s adequately ‘recharged. | ain rs 
_ Problems of inadequate transmissibility of aquifers have developed in every 
part of the country. - Pumping from a well or closely spaced group of wells — 
we 
creates a cone of depression in the water table or pressure surface of an i 
aquifer. Generally, if withdrawal of water continues at a constant rate, the : 

"cone ¢ expands and the pumped water draws water from a progressively increas- a 
ing area. . The v water level in the well continues to drop, but at a decreasing i 
‘rate, until me cone has expanded to reach either an area of natural recharge “- 

: ord discharge. © At such time as the amount withdrawn is balanced by increased _ 
4 ‘movement from the recharge a area or decreased natural discharge, the 5 pro- 7 
gressive lowering of the water table will cease due to stabilization of rates. = 
_ The ‘solution | to problems of apparent | or local shortages of water is to a 
fect a balance between the rate of draft and the rate of replenishment, either 
by reducing the draft or increasing the replenishment, or both. The principle 
— the same as that involved in the elimination of overdraft from reservoirs — 
where long accumulated storage is being progressively depleted. Other cor- 
_ rective measures that are effective in many areas include development ll 
wells in previously untapped portions of the aquifer or redistribution of wells 
to draw from a more extensive part of the aquifer. Such measures tend to | 


reduce concentrated draft upon a small portion of the aquifer. 


_ This brief resume of the general nature of problems connected with ground © 
2 water development will, it is hoped, serve to indicate the intricate nature of Bi 3 


considerations affecting a given situation. These considerations are 
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Discussion by. Fred W. Blaisdell and Philip W. Manson 


FRED W. BLAISDELL, 1 ASCE ont PHILIP Ww. MANSON.2 —The 


= done an excellent job of analyzing and presenting the mass of data an 7 


‘The writers conducted tests on Square- edged pipe junctions between 


stream pipes were of the same size and this discussion will be restricted to 
in-line pipes of the same size. The writers’ laterals covered a wider range 


are different in form. In the writers’ work, ‘the upstream and 


_ Of sizes than did the authors’. Whereas the authors tested only a 90° junction, 
the writers’ junction angles varied from 15° to 165° in 15° increments. Only 
“the 90° junction will be compared with the authors’ results. In the > 
techniques were used by both experimenters. 
spite of the large number of tests performed by each experimenter, a _ 
dir 


ect comparison of the results is possible only for the one case in which all © 


three pipes are of the same size. The in-line loss coefficient agrees very 
well with data obtained by the writers and other experimenters for the square — 
junction, and less well for the rectangular junction box when a large portion — 
of the total flow enters from the lateral. Similar good agreement between the 


‘lateral loss coefficient except when all the flow is from the lateral. 
7 Rke It is reasonable to expect better agreement between the authors’ and a 
writers’ results for the square junction than for ' the rectangular junction be- 


by ‘the writers. 
a the agreement of the two sets of data | is remarkable. 
Simple reasoning would indicate that Eq. (4) whith applies so 
well to the rectangular junction, should apply to square and round junctions _ 


there a an explanation why different formulas were found 


a. Proc. Pape Paper 2057, 1959, by M. Horace W. Wood, 
Ernest T. Smerdon and HerbertG. Bossy, 


8. De 


_ Project Supervisor, ‘Agri. Research Service, U. S. Dept. of | Agriculture, 


St. Anthony Falls Hydr. Lab., Minneapolis, Minn. 
Prot. Agri. Univ. of Minn. , St. Paul, Minn. 
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